M D /"O 6 

/ 

i CORRELATION OF TEXTURES AND 

MICROSTRUCTURES IN A SERIES OF 

Ni-Co ALLOYS 


By 

Krishnendu Mukherjee 







DEPARTMI|IT OF MATERIALS AND METALLURGICAL ENGINEERING 


Indiar 


institute 


of Technology 

APRIL, 2002 


Kanpur 


CORRELATION OF TEXTURES AND 
MICROSTRUCTURES IN A SERIES OF 

Ni-Co ALLOYS 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 


MASTER OF TECHNOLOGY 


By 

KRISHNENDU MUKHERJEE 


To the 


DEPARTMENT OF MATERIALS AND METALLURGICAL ENGINEERING 

Indian Institute of Technoiogy, Kanpur 


APRIL, 2002 



/ 

""4 FEB 2M 

^^rfccT 3Eio A *mymm iMte <Mw 



CERTIFICATE 


This is to certify that the present work, entitled “CORRELATION OF 
TEXTURES AND MICROSTRUCTURES IN A SERIES OF Ni-Co ALLOYS” by 
Mr. Krishnendu Mukheijee has been carried out under our supervision and to the best of 
our knowledge it has not been submitted elsewhere for a degree. 



Dept, of Materials & Metallurgical Engg. 
Indian Institute of Technology 
Kanpur 


Date: 


0 ^ 



Dr. S. Sangal 

Professor 

Dept, of Materials & Metallurgical Engg. 
Indian Institute of Technology 
Kanpur 


Date: 



ACKNOWLEDGEMENT 


I would like to convey my deepest sense of indebtedness and gratitude to my thesis 
supervisors Prof. R. K. Ray and Prof. S. Sangal for their inspiring guidance in every 
phase of this work. 

I sincerely acknowledge the help and support rendered by Dr. Mungole to allow me to 
use the Image Analyzer of the Metallography Lab. 

I would like to convey my sincere thanks to Mr. V. Kumar for extending all possible and 
sincere help through the course of this study and I feel very fortunate to have a laboratory 
assistant like him. 

Finally I would like to mention those integral part of my IITK life, my friends, who have 
made every moment of stay here at IITK to rejoice in future. The list is endless, still, the 
support, cooperation and help rendered by my laboratory partners Suhas, Pinaki, Arijit, 
Deepti, Lt. Sanjay Chhabra, S. Sankaran, John, Subhashis, Kausik, and Pushpender 
Chhillar are gratefully acknowledged. Also I am thankful to Debu, Tapas Laha, Sudin, 
Natida, Sona, Amlan, Parthada, Sagnik, Santanu, Dulal, and Projitda for helping me in 
various ways. 


Krishnendu Mukherjee 




CONTENTS 


Page No. 

LIST OF FIGURES 
LIST OF TABLES 

ABSTRACT 

1. INTRODUCTION 1 

2. LITERATURE REVIEW 4 

2.1 THE NICKEL-COBALT SYSTEM 4 

2.2 THE COLD WORKED STATE 4 

2.2.1 MICROSTRUCTURAL FEATURES OF THE COLD WORKED 

STATE 6 

2.2.2 THE STORED ENERGY PRODUCED BY COLD WORKING 6 

2.3 RECOVERY 6 

2.3.1 RECOVERY MECHANISM 6 

2.4 PRIMARY RECRYSTALLISATION 7 

2.4.1 CLASSICAL NUCLEATION MODEL 7 

2.4.2 SUBGRAIN GROWTH MODEL 7 

2.4.3 STRAIN INDUCED BOUNDARY MIGRATION 7 

2.5 GRAIN BOUNDARY MIGRATION 8 

2.6 CONTINUOUS GRAIN GROWTH 9 

2.7 MEASUREMENT OF TEXTURE 9 

2.8 ROLLING TEXTURE OF Ni-Co ALLOYS 14 

2.9 RECRYSTALLISATION TEXTURE 18 

2.10 ESTIMATION OF GRAIN BOUNDARY ENERGY DISTRIBUTION 22 

2.1 1 RELATION BETWEEN PDA AND TDA 24 



3. MATERIAL, MICROSCOPY AND IMAGE ANALYSIS 


,32 


3.1 MATERIAL AND INITIAL TREATMENT 32 

3.2 COLD ROLLING 32 

3.3 ANNEALING 32 

3.4 TEXTURE MEASUREMENT 33 

3.5 OPTICAL MICROSCOPY 33 

3.6 IMAGE PROCESSING ALGORITHM 33 

3.6.1 BINARIZATION/ THRESHOLDING 34 

3.6.2 DILATION FILTERING 34 

3.6.3 EROSION FILTERING 34 

3.6.4 THINNING MORPHOLOGICAL FILTER 34 

3.6.5 BITMAPS 40 

3.7 TRIPLE POINT DETECTION 40 

3.8 EVALUATION OF PLANE DIHEDRAL ANGLES (PDA) 43 

3.9 TRANSFORMATION OF PDA DISTRIBUTION TO TDA 

DISTRIBUTION 46 

3.10 RELATIVE GRAIN BOUNDARY ENERGIES 46 

4. RESULTS AND ANALYSIS 47 

5. DISCUSSION 122 

6. CONCLUSIONS 125 

REFERENCES 126 


APPENDIX 



LIST OF FIGURES 


Figure 

No. 


2.1 


2.2 














Title of the figures 


Nickel-Cobalt equilibrium diagram 

The rolling plane and rolling direction in a metal sheet 

subjected to rolling. 

Schematic illustration of pole figure construction 


(a): Specimen frame and crystallite frame. 

(b): The rotations for superimposing crystallite frame on 
reference frame. 

(c): Orientation space. 


Effect of Composition on Stacking Fault Energy 


{111} Pole figures showing the rolling textures of: (a) pure Ni; 
(b) Ni-10% Co; (c) Ni-20% Co; (d) Ni-30% Co; (e) Ni-40% 
Co; and (f) Ni-60% Co. 


Plot Showing the Variation of M, (Volume Fraction) with % Co 
in Ni-Co alloys. 


Texture of pure Copper (a) Cu type rolling texture, (b) Cube 
type annealing texture. 


Schematic representation of growth of Cube Texture region 
during annealing. 


Sectioning of a triple edge by a plane of polish randomly 
oriented at (0, (j)). ai, a 2 , and as are the true dihedral angles and 
Pi, P 2 and Ps are the corresponding plane dihedral angles; yi, ya 
and ys are the energy per unit area of the grain boundary planes. 


An optical microstructure of a polycrystalline material, pi, p 2 
and p 3 are the plane dihedral angles at a triple point. 


A polycrystalline structure (X) consisting of randomly oriented 
triple edges being sectioned by the plane of polish (T 2 ). 


A random section through two grain boundary planes meeting 
at an angle of a. The normal to the plane of section makes an 
angle 0 with the z axis and the projection of the normal on the 
XY plane makes an angle (j) with the X axis. The line of 
intersection of grain boundaries on plane of polish makes an 

angle p. 


(a) A typical grayscale image of a polycrystalline 
microstructure, (b) Bimodal gray level histogram of the image 

(a). 


Transformed binary image (after applying thresholding 
algorithm) showing polishing/etching artefacts. 


Result of applying dilation filter on Figure 3.2. The 
disappearance of peninsulas and etching artefacts is noticeable. 


Image of microstructure obtained after erosion filter was 
applied on the image shown in Figure 3.3. It is noticeable here 
that the grain boundary has become continuous. 












































3.5 Image of the resultant microstructure after application of 39 

thinning filter on the image shown in Figure 3.4. The grain 
b ounda ries are now o nl y on e pii;^[thiclc 


3.6 Bitmap of the marked rectangular region in Figure 3.5. ‘0’ 41 

refers to grain boundary pixels and ‘ 1 ’ refers to matrix pixels. 


A partial configurations of valid triple points. Remaining 42 

configurations can be generated rotating the matrices through 

7r/2, n and 3m/2. 


3.8 Illustration of pseudo triple points. ^ 


3.9 Tangents (thick lines) drawn to the fitted grain boundary 45 

segments (thin lines) at a triple point ‘O’. ‘A’ and ‘B’ are 
arbitr ar y points on the segments. 


3.10 Schematic of cumulative distribution for true dihedral angles: 45 

(a) amiis . (h) ainm- Pi and P 2 arc randomly generated number. xi 
and X 2 are corresponding Umax and amm- 


4l I Micrographs of different Ni-Co alloys annealed for 1 5 minutes. 48 


.2 Micrographs of different Ni-Co alloys annealed for 30 minutes. 


3 Micrographs of different Ni-Co alloys annealed for 1 hour. 


Micrographs of different Ni-Co alloys annealed for 3 hours. ^ 


5 Micrographs of different Ni-Co alloys annealed for 10 hours. ^ 


.6 r Micrographs of different Ni-Co alloys annealed for 20 hours. 53 


Micrographs of different Ni-Co alloys annealed for 50 hours. 54 


4.8 Micrographs showing processed structure of different Ni-Co 55 

alloys annealed for 1 hour. 


Micrographs showing processed structure of different Ni-Co 56 
alloy annealed for 50 hours. 


4. 1 0 Frequency distribution of relative energy ratio yl /yS of different 58 
Ni-Co alloys annealed for 15 minutes. 


4. 1 1 Frequency distribution of relative energy ratio y l/y3 of different 59 

Ni-Co alloys annealed for 30 minutes. 


4.12 Frequency distribution of relative energy ratio yl/yl of different 60 

Ni-Co alloys annealed for 1 hour. 


4.13 Frequency distribution of relative energy ratio yl/y3 of different 61 

Ni-Co alloys annealed for 3 hours. 


4.14 Frequency distribution of relative energy ratio yl/y3 of different 62 

alloys annealed for 10 hours. 


4. 1 5 Frequency distribution of yl/y3 energy ratio of different alloys 63 

annealed for 20 hours. 


4. 1 6 Frequency distribution of relative energy ratio yl/y3 of different 64 

Ni-Co alloys annealed for 50 hours. 


4. 1 7 Frequency distribution of relative energy ratio y2/y3 of different 65 
Ni-Co alloys annealed for 15 minutes. 


4. 1 8 Frequency distribution of relative energy ratio y2/y3 of different 66 

Ni-Co alloys annealed for 30 minutes. 


4. 19 Frequency distribution of relative energy ratio y2/y3 of different 

Ni-Co alloys annealed for 1 hour. 


4.20 I Frequency distribution of relative energy ratio y2/y3 of different I 68 
Ni-Co alloys annealed for 3 hours. 




















































4.21 

Frequency distribution of relative energy ratio y2/y 3 of different 
Ni-Co alloys annealed for 10 hours. 

69 

4.22 

Frequency distribution of relative energy ratio yllyi of different 
Ni-Co alloys annealed for 20 hours. 

70 

4.23 

Frequency distribution of relative energy ratio yllyl of different 
Ni-Co alloys annealed for 50 hours. 

71 

4.24 

Frequency distribution of relative grain boundary energy ratio 
y1/y 3 for Ni-10% Co samples subjected to different annealing 

times. 

72 

4.25 

Frequency distribution of relative grain boundary energy ratio 
Y1 /y 3 for Ni-20% Co samples subjected to different annealing 

times. 

73 

4.26 

Frequency distribution of relative grain boundary energy ratio 
y1/y 3 for Ni-30% Co samples subjected to different annealing 

times. 

74 

4.27 

Frequency distribution of relative grain boundary energy ratio 
y1/y 3 for Ni-40% Co samples subjected to different annealing 

times. 

75 

4.28 

Frequency distribution of relative grain boundary energy ratio 
y1/y 3 for Ni-60% Co samples subjected to different annealing 

times. 

76 

4.29 

Frequency distribution of relative grain boundary energy ratio 
Y2/y3 for Ni-10% Co samples subjected to different annealing 

times. 

77 

4.30 

Frequency distribution of relative grain boundary energy ratio 
y2/y 3 for Ni-20% Co samples subjected to different annealing 

times. 

78 

4.31 

1 

Frequency distribution of relative grain boundary energy ratio 
y2/y 3 for Ni-30% Co samples subjected to different annealing 

times. 

79 

i 

4.32 

Frequency distribution of relative grain boundary energy ratio 
y2/y 3 for Ni-40% Co samples subjected to different annealing 

times. 

80 

4.33 

Frequency distribution of relative grain boundary energy ratio 
y2/y 3 for Ni-60% Co samples subjected to different annealing 

times. 

81 

4.34 

The variation of frequencies of different classes of yMyZ as a 
function of annealing time for Ni-10% Co samples. 

83 

4.35 

The variation of frequencies of different classes of y1/y 3 as a 
function of annealing time for Ni-20% Co samples. 

: 84 

4.36 

The variation of frequencies of different classes of yMylt as a 
function of annealing time for Ni-30% Co samples. 

85 

4.37 

The variation of frequencies of different classes of y1/y 3 as a 
function of annealing time for Ni-40% Co samples. 

86 

4.38 

The variation of frequencies of different classes of yMyT) as a 
function of annealing time for Ni-60% Co samples. 

87 

4.39 

The variation of frequencies of different classes o{y2ly'i as a 
function of annealing time for Ni-10% Co samples. 

88 


































4.40 


89 


4.41 


4.42 


4.43 


4.44 


4.45 


4.46 


4.47 


4.48 


4.49 


4.50 


4.51 


4.52 


4.53 


4.54 


4.55 


4.56 


4.57 


4.58 


4.59 


The variation of frequencies of different classes of Y2/y3 as a 

function of annealing time for Ni-20% Co samples. 

The variation of frequencies of different classes of Y2/y3 as a 

function of annealing time for Ni-30% Co samples. 

The variation of frequencies of different classes of y2/y3 as a 
function of annealing time for Ni-40% Co samples. 

The variation of frequencies of different classes of y2/y3 as a 

function of annealing time for Ni-60% Co samples. 

The variation of mean of relative energy ratio yl/y3 of different 

alloys as a function of annealing time. 

The variation of mean of relative energy ratio y2/y3 of different 

alloys as a function of annealing time. 

The variation of co-cfficicnt of variation of relative energy ratio 
yl/y3 of different alloys as a function of annealing time. 

The variation of co-efficient of variation of relative energy ratio 
y2/y3 of different alloys as a function of annealing time. 
Variation of intensity of cube component as a function of 

annealing time for different samples. 

Variation of volume fraction of cube component as a function 

of annealing time for different samples. 

Variation of intensity of cube component as a function of Co 

content and subjected to different annealing times. 

Variation of volume fraction of cube component as a function 
of Co content and subjected to different annealing times. 
Variation of volume fraction of cube and its twin component as 

a function of annealing time for different samples. 

Variation of volume fraction of cube and its twin component as 
a function of Co content and subjected to different annealing 

times. 

Variation of intensity of BsT component {437}<4 11 7> as a 

function of annealing time for different samples. 

Variation of volume fraction of BsT component {437}<4 1 1 7> 
as a function of annealing time for different samples. 
Variation of intensities of different texture components in Ni- 

40% Co alloy as a fimetion of annealing time. 

Variation of volume fractions of different texture components 

in Ni-4Q% Co alloy as a function of annealing time. 

Variation of intensities of different texture components in Ni- 

60% Co alloy as a function of annealing time. 

Variation of volume fractions of different texture components 
in Ni-60% Co alloy as a function of annealing time. 


90 


91 


92 


93 


94 


96 


97 


108 


109 


111 


112 


113 


114 


115 


116 


118 


119 


120 


121 










































LIST OF TABLES 










ABSTRACT 


An attempt has been made to correlate crystallographic texture with microstructure in 
a series of Ni-Co alloys in the recrystallised condition and during the grain growth 
stage. For this purpose Ni-10%Co, Ni-20%Co, Ni-30%Co, Ni-40% Co, and Ni- 
60%Co (all in wt.%) alloys were cold rolled 95% at room temperature and then 
subsequently recrystallisation annealed at 800“C for various lengths of time varying 
between 15 minutes to 50 hours. Crystallographic texture in the form of pole figures 
and ODFs, were measured for the different heat-treated samples. Optical microscopy 
was carried out and an automated methodology, based on image processing and 
analysis, has been developed in order to evaluate the plane dihedral angles and then 
the true dihedral angles in the microstructures. The relative grain boundary energy 
ratios were evaluated from the true dihedral angles. 

The recrystallisation textures of Ni-10% Co to Ni-30% Co alloys are very similar. 
These are comprised of mainly the cube, its twin components, and the BsT 
component.In Ni-40% Co alloy the important recrystallisation texture components are 
the BsR component, the {279}<7 11 7> component, the retained S component, and 
the cube component.The recrystallisation texture of Ni-60% Co alloy consists mainly 
of rotated Goss component and the component {011}<114>. The frequency 
distribution of the grain boundary energy ratios correspond well with the volume 
fraction of the texture components in all the five alloys after 20 hours of annealing 
treatment. 



CHAPTER 1 

INTRODUCTION 


The properties of a polycrystalline material are functions of both microstructure and 
grain orientation. By varying the processing parameters of a material we can get 
desired microstructurc as well as grain orientation. If a correlation between the two 
can be found out we can get both desired microstructurc and desired grain orientation 
together. There are different aspects of microstructure e.g. grain size, grain boundary 
roughness, dihedral angle. Amongst them dihedral angle is an important parameter as 
it gives an idea of the relative grain boundary energy. 

The common materials of industrial practice are polycrystalline aggregates in which 
each of the individual grains may have an orientation that differs from those of its 
neighbors. In natural practice however, it is quite unusual for the grains in such 
materials to have a random distribution of orientations. These non-random 
distributions are called preferred orientations or texture. 

Texture in material is developed by mainly four ways: 

1 . Solidification texture: It is the texture developed during the solidification of liquid 
material. 

2. Deformation texture: A material which has undergone a severe amount of 
deformation as in rolling will develop a preferred orientation, or texture, in which 
certain crystallographic planes tend to orient themselves in a preferred manner 
with respect to the direction of maximum strain. In case of a single crystal the slip 
planes tend to rotate parallel to the axis of principal strain. The same situation 
exists in a polycrystalline aggregate, but the complex interactions between the 
multiple slip system makes analysis of the polycrystalline situation much more 
difficult. Since the individual grains in a polyciystalline aggregate cannot rotate 
freely, the lattice bending and fragmentation will occur [Ref. 1]. 

3. Annealing texture: T he recrystallisation of a cold-worked material generally 
produces a preferred orientation which is different from that existing in the 
deformed material. This is called an annealing texture, or recrystallisation texture. 
An outstanding example is the cube texture in copper, where the {100} plane lies 
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parallel to the rolling plane with a <001> direction parallel to the direction of 
rolling. 

4. Transformation texture: It is the texture developed during the change of phase by 
heat treatment. 

Deformation texture is an internal property of the material itself It depends on the 
deformation mode of the material and the deformation mode depends on the stacking 
fault energy (SFE) of the material, which is a material property. 

In FCC materials usually three different types of deformation textures are observed: 

1 . Brass or alloy type textures in low SFE materials. 

2. Copper or pure metal type texture in medium SFE materials. 

3. Aluminium type texture in very high SFE materials. 

While in the brass type there is only the brass component {110}<112> but in the 
copper type the Cu {112}<111>, S {123}<634> and brass components are nearly 
equally strong. The S component predominates in the aluminium type texture. 

Pure nickel has high SFE (~ 250 mJ m‘^). Addition of cobalt to nickel reduces the SFE 
of the alloy drastically. So, in the Ni-Co binary (FCC solid solution) a wide range of 
alloys with different values of SFE can be developed just by changing the 
composition (i.e. the cobalt content). Change of composition directly affects the 
deformation mode and thereby deformation texture. This is also going to be reflected 
in the primary recrystallisation textures in these alloys. The textural changes are found 
to persist even after primary recrystallisation, i.e. during the continuous grain growth 
stage. Since grain growth essentially consists of the growth of crystallographically 
favorably oriented grains at the expense of others, grain growth should be intimately 
related to the textural changes. 

When a high or medium SFE FCC material is annealed after heavy cold working, a 
strong cube texture {100}<001> forms. This is a familiar phenomenon for the pure 
FCC metals like Cu, Ni etc. Aluminium is an exception, where the rolling texture is 
dominated by S type texture, and on recrystallisation along with the cube texture, R 

texture (retained rolling texture) {123}<412> is obtained which is very close to S 
type texture. On the other hand, on annealing the Bs type rolling texture yields a 
series of texture components. 

The properties of polycrystalline materials depend a lot on the grain boundary energy. 
The grain boundary energy state influences the strengthening properties of the 
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polycrystals significantly. True dihedral angle i.e., angle between grain boundary 
planes meeting at a triple edge in three-dimensional structure of polycrystals is 
indicator of the energy state of the grain boundaries. True dihedral angle distribution 
as a niicrostructural parameter has not been generally used because of the 
experimental difficulty in its measurement. Flowever a recently developed technique 
[Ref, 2 \ makes it possible for the measurement of true dihedral angles rather simple. 
This technique measures the true dihedral angles from the plane dihedral angles, 
which can be determined comparatively easily. Plane dihedral angles are measured 
from the microstructure with the help of suitable computational technique involving 
digital data processing. The application of digital data processing has several 
advantages e.g., statistically adequate number of data can be processed in reasonable 
time, accuracy level is increased. From the true dihedral angle distribution relative 
energies of grain boundaries can be determined. Grain boundary energy ratios can be 
easily determined and this can be correlated with the textural aspects of the material. 


CHAPTER 2 

LITERATURE REVIEW 


2.1: THE NICKEL-COBALT SYSTEM 

In the periodic classification of elements, Nickel and Cobalt both occur in the first 
long period (period 4) of the transition group VIII [Ref. 3]. The basic physical 
properties of Nickel and Cobalt are depicted in Table 2.1. 

TABLE 2.1: PHYSICAL PROPERTIES OF NICKEL AND COBALT 



Nickel 

Cobalt 

Symbol 

Ni 

- , 

Co 

Atomic Number 

28 

• . ' 27 • • 

Atomic Weight 

58.71 

58.93 

Crystal Structure 

FCC 

HCP 

Density 

8908 kg m’^ 

8900 kg m'^ 

Melting Point 

1453'C 

1495'C 

Boiling Point 

2730“C 

2927°C 

Electronic Structure 


[Ar]3d'^4s^ 

Lattice Parameter 

3.5238 A' 

a=2.5071 A" 

b=2.5071 A' 

c=4.0695 A* 


• [Ar]=ls^2s^2p®3s23p^ 


Figure 2.1 shows the diagram of the nickel cobalt system [Ref. 4J. The effects of 
cobalt additions to nickel are to increase the Curie point and the magnetic saturation, 
and to slightly improve the mechanical properties. Nickel-cobalt alloys containing 
upto about 30% cobalt are formed by electrodeposition from mixed-salt electrolytes, 
and they have an exceptionally high hardness. 

2.2: THE COLD WORKED STATE 

The cold worked state may be taken as any strained or damaged condition of a 
crystalline material brought about by processes such as plastic deformation, particle 
bombardment, quenching from a high temperature or phase transformations. [Ref. 5] 
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The principal entities responsible for the cold worked state are point defects such as 
vacancies and interstitials as well as linear and planar defects like dislocations and 
stacking faults etc. If the deformation occurs at a sufficiently low temperature, many 
of these defects will be retained, leading to a thermodynamically metastable state. 
2.2.1: MICROSTRUCTURAL FEATURES OF THE COLD WORKED STATE 
The regions inside the grains of deformed material are broken up into units of slightly 
misoriented “crystal domains”, which are called cells. The cell interiors are relatively 
free from dislocations, whereas the walls are regions of high dislocation density. The 
dislocation arrangements within the cell walls are complex and are often associated 
with small loops. The misorientations across the cell walls are found to be of the order 
of a degree. An increased amount of deformation results into the increase in 
misorientation even upto 10°. 

2.2.2: THE STORED ENERGY PRODUCED BY COLD WORKING 
A small percentage (usually from 1% to 10%) of the energy expended in plastically 
deforming a material remains stored in the material causing an increase in the internal 
energy. The magnitude of the stored energy produced by cold working is affected by 
variables such as purity, method of deformation, temperature, grain size and 
composition. The energy of plastic deformation is stored in the following form: 

1. Elastic strain energy 

2. Energy of dislocation 

3. Energy of point defects 

4. Energy associated with stacking faults and twins 

5. Energy associated with the destruction of ordered region 
2.3: RECOVERY 

Recovery includes the migration of vacancies and of dislocations, leading to 
annihilation and rearrangement of certain portions of the grain. Recovery can be 
considered to be any modification of properties during annealing, which occurs before 
the appearance of the new strain free recrystallised grains. 

2.3.1; RECOVERY MECHANISM 

A cold worked metal has a cell structure. As annealing starts, the tangled dislocations 
in the cell wall are supposed to rearrange themselves. Some of the dislocations within 


the cells are attracted to cell walls. As a result the dislocation density in the cell 
interior decreases gradually. The cell walls become more clearly defined and 
eventually form sub-grains of about the same size as the initial cells. At the later stage 
it starts increasing gradually. 

2.4: PRIMARY RECRYSTALLISATION 

Recrystallisation of a deformed metal involves the nucleation and growth of new 
strain free grains, which may be widely different in orientation from the surrounding 
matrix. The new grains grow at the expense of the cold worked matrix by high angle 
boundary migration. The three principal models of nucleation are: 

1 . Classical nucleation model, 

2. Sub-grain growth model, 

3. Strain induced boundary migration model. 

2.4.1: CLASSICAL NUCLEATION MODEL 

This is based on the assumption that the nucleus, which is stable and capable of 
growth, is formed as a result of thermal fluctuations in the sense of the classical 
nucleation theory of the phase transformation. This model requires the existence of an 
incubation period. It also predicts the preferential nucleation of recrystallised grains in 
the most severely deformed regions. But the main limitation of this theory is that it 
predicts the new grains should always be in close orientation with the adjacent cold 
worked matrix. Another drawback is that there is no consideration for purity of the 
material. However, in practical situations this theory is quite applicable. 

2.4.2: SUBGRAIN GROWTH MODEL 

Small strain free cells can be formed in a deformed matrix by dislocation climb and 
rearrangement. Once a sufficiently large dislocation free cell or sub-grain has formed 
it can grow into its neighbourhood. The misorientation of a sub-grain relative to its 
immediate surrounding is more important than its size in controlling its growth. 

2.4.3: STRAIN INDUCED BOUNDARY MIGRATION 

This is basically a model of growth since it assumes that the unstable cold worked 
state already has large angle boundaries of macroscopic or microscopic dimensions 
and that during recrystallisation these boundaries simply bow out leaving strain free 
materials behind. Sometimes recrystallisation does not involve any nucleation but 
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only growth of strain free regions brought about by grain boundary migration. This 
process which has been termed as strain induced boundary migration, consists of the 
irregular movement of an existing boundary with an irregular jagged front. According 
to this model there is no incubation period, since no nucleus of new orientation has to 
be formed. It is observed that the victim grain is the one which is more severely work 
hardened, as the stored energy of this grain provides the driving force for the entire 
process. The grain with larger cell size will have the lower energy, so that 
reciystallisation occurs by the movement of a portion of the boundary into the high 
energy grain having a smaller cell size. 

2.5: GRAIN BOUNDARY MIGRATION 

Grain boundary migration plays an essential role in recrystallisation and subsequent 
grain growth. There is a marked orientation dependence of grain boundary migration. 
The highest migration rate in FCC metals has been obtained when the growing grain 
is crystallographically related to the matrix by a 30° - 40° rotation about a common 
<1 1 i> axis [Ref. 6). 

The anisotropy of the boundary migration rate can be interpreted as due to the 
orientation dependence of either the driving force or the boundary mobility. In highly 
deformed materials, which show a complex dislocation structure, the anisotropy of 
driving force can be ruled out as a possibility and hence an orientation dependence of 
the boundary mobility is quite expected. 

Reciystallisation was found to be easier In foils made from the transverse and 
longitudinal sections than from ones parallel to the rolling plane [Ref. 7]. This is due 
to the fact that the ceils produced by deformation are pancake shaped with the cell 
boundaries lying parallel to the rolling plane. Foils, which were taken parallel to this 
plane, may not contain a sufficient number of ceils through the thickness to allow 
formation of a recrystallisation nucleus with a high angle boundary. In such foils 
growth may also be inhibited. 

The kinetics of secondary recrystallisation is affected by 

1 . The value and behavior of pinning force, 

2. The magnitude of the driving force for grain growth [Ref. 8]. 
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The former depends on the volume fraction, size, and the nature of the fine 
precipitates. The later is determined by. initial grain size distribution and curvature of 
the moving boundaries. The two dimensional growth slows down the secondary 
recrystallisation when the initial average grain size exceeds one-tenth of the sheet 
thickness. 

2.6: CONTINUOUS GRAIN GROWTH 

Continuous grain growth is a recrystallisation process, which proceeds without 
nucleation. Textures formed during this process are thus exclusively due to the 
orientation dependence of the grain boundary velocity which is the product of the 
driving force and grain boundary mobility. The mobility of a boundary depends 
among others, on the orientation difference Ago of the contiguous grains. The driving 
force in the grain coarsening is constituted by the deviation from equilibrium between 
the forces acting on the grain boundaries. These forces also depend on the orientation 
difference of the contiguous grains. In order to reduce deviation from the equilibrium 
state, the system lowers its free energy by enlarging the grains, thereby reducing the 
surface free energy. So, this is essentially a process, where a grain of a favoured 
orientation consumes other grains of unfavorable orientations. 

2.7: MEASUREMENT OF TEXTURE 

Rolling plane of the material and rolling direction (crystallographic planes and 
direction) are needed for the determination of texture. For example, when a sheet is 
rolled in such a way that the (hkl) planes of most of the grains are parallel to the 
rolling plane, and the [uvw] direction of most of the grains is parallel to the rolling 
direction, then texture of the rolled sheet is represented as {hkl}<uvw>, as shown in 
Figure 2.2. This is displayed by a pole figure, which is nothing but the stereographic 
projection on the (hkl) plane, taking (uvw) pole as the north pole. The experimental 
procedure of the pole figure construction is briefly 

shown in Figure 2.3. 'I'he material is first subjected to X-ray diffraction, and the 
intensity vs. angle data is recorded. Then these data is transformed from a linear scale 
to a polar plot (pole figure). Several contours shown on the pole figure account for the 
different orientations of different densities. 
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ClOO> poles for the orientation {I10)<00l> 










There are few limitations of pole figure. Although exhibiting all the texture 
components of a material, in a pole figure, the pole of particular texture is not 
pinpointed, rather it is scattered throughout a zone and it becomes difficult to 
distinguish them perfectly. The resolution of pole figure is poor. Also, it gives only 
some qualitative idea about the strength of a particular texture component. 

One-way of removing this difficulty is to use the crystallize orientation distribution 
function (CODF or ODF^, which essentially describes the frequency of occurrence of 
particular orientations in three dimensional (Eular) orientation space. It is easily 
understood that a frame of three mutually perpendicular axes can be superimposed on 
another frame of three mutually perpendicular axes. The same principle is followed 
here, the specimen frame as shown in Figure 2.4(a) is superimposed on the 
crystallographic frame by means of three successive rotations (Figure 2.4(b)). 

(a) A first rotation around ND transforms TD and RD into the new 
directions TDi and respectively. Oi has to have such a value that RD| will be 
perpendicular to the plane formed by ND and [001]. 

(b) A second rotation <I) around the new direction RDj with O having such a 
value that ND is transformed into [001] and TD into TDj. 

(c) The third rotation <I >2 with O 2 having such a value that RDj is transformed 
into [100] and TD 2 into [010]. 

So, Oi,C>,02; the set of these three angles completely represent an orientation of a 
crystal or grain. Immediately for every set of {hkl}<uvw> values there is exactly one 
set of Oi O O 2 values. Now depending on this concept, a three dimensional space has 
been evaluated where three mutually perpendicular axes represent the values of 
01 , 0 , 02 ; from 0° to 90°. This is known as orientation space (Figure 2.4(c)), and 
widely used for universal texture notation [Ref. 9]. 

The ODF technique is not more correct process than the pole figure technique is, but 
of course, the resolution is much higher in this system, and also it allows the 
quantitative measurement of the strength of a particular texture component in volume 
fraction. 




RD 

SperiM» O f™** 



Figure 2.4(a): Specimen frame and crystallite frame. 

F^re 2.4(b): The rotations for superimposing crystallite frame on reference frame. 
Fi^re 2.4(c): Orientation space. 
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2.8: ROLLING TEXTURE OF Ni-Co ALLOYS 

For higher stacking fault energy, the stacking fault width is very low, and the stacking 
faults are easily constricted. In that situation, the cross-slip becomes easier and the 
material undergoes slip deformation rather than twinning. When the deformation is 
governed by slip mode, the deformation texture is generally copper type { 1 12}<1 1 1>. 
Example, pure copper, pure nickel or nickel with low amount of cobalt, Ni-Fe alloy 
etc. Along with the copper components of texture other components like Bs 
component {110}<112>, S component {123}<634>, Goss component {110}<001> 
^ etc. are present in very small percentage. 

For lower SFE, the stacking fault width is much greater, which makes difficult to 
create any constriction within the stacking fault region and the cross-slip becomes 
difficult. As a result of it, the critical resolved shear stress (CRSS) for twin becomes 
lower than that of slip and the plastic deformation is essentially governed by the 
twinning process. 

In Ni-Co system, the SFE significantly decreases on increasing the Co percentage 
(Figure 2.5). Accordingly it is observed that the Cu component of rolling texture 
gradually decreases from pure nickel to Ni-30% Co (Figure 2.6). Ni-40 Co shows a 
transition from Cu type to Bs type rolling texture, and Ni-60 Co alloy predominantly 
exhibits Bs type texture [Ref. 10]. 

This is also observed that Cu type texture decreases with Increasing the amount of 
deformation. In case of Ni-40% Co alloy, the intensity of Cu type texture remarkably 
decreased from 40% cold rolling to 95% cold rolling. 

Wassermann has shown that if mechanical twinning is an available mode of 
deformation in addition to normal slip, the Cu component of the rolling texture 
{1 12}<1 1 1> will be transformed to {552}<115> by twinning. Further deformation of 
the twin orientation will then take place on the favorable slip systems which will 
rotate this twin orientation into the Goss position 
{0I1}<100> and finally into the Bs position {0I1}<112> by further rotation. In 
contrast, the major component of the alloy type texture will not be affected by 
twinning and retain its orientation during deformation. Figure 2.7 clearly indicates 
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Figure 2.5: Effect of Composition on Stacking Fault Energy 
(y = Stacking Fault Energy, G = Shear Modulus, b = Burger’s Vector) 
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Figure 2.7; Plot Showing the Variation of Mj (Volume Fraction) with % Co in Ni-Co 
alloys 
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the phenomenon of dominating Cu type texture in lower Co content alloys and that Bs 
texture in higher Co content alloys. 

Temperature plays an important role for determining the rolling texture. The slip 
systems are activated in larger extent at high temperature. So, at high temperature, slip 
deformation is more probable than twinning. Exactly, the reverse effect is expected at 
lower temperatures. 

On annealing a deformed material, the rolling texture completely or partially changes 
to a new texture, which is the recrystallisation or annealing texture. It is also 
observed, that there is a one to one correspondence between the rolling texture and the 
annealing texture. So it can be inferred that the ultimate annealing texture of the 
material is dependent on the following factors: 

1 . Deformation mode, 

2. Degree of deformation, 

3. Temperature at which rolling is carried out, 

4. Stacking fault energy, 

5. Critical resolved shear stress for slip and twinning. 

2.9: RECRYSTALLISATION TEXTURE 

When the cold worked material is annealed, the Cu type texture essentially changes to 
cube texture {100}<001>. This is a familiar phenomenon for the pure FCC metals 
like Cu, Ni etc. The reverse effect is observed in Bs type rolling texture. On 
recrystallisation it does not yield Cu type texture. The Bs texture transforms to 
{332}<1 1 3> type texture. 

The formation of cube texture is not a desirable one so far as the material property is 
concerned. The formation of cube texture is a very fast process, and just after the 
starting of recrystallisation the complete material immediately becomes highly 
anisotropic. But the reason for the dominating characteristic of cube texture and its 
actual origin is still not very clearly known. The pole figures of 95% cold rolled pure 
copper and the recrystallised texture of the same have been juxtaposed in the Figure 
2.8 from which it is clearly observed that in the cold deformed alloy, a high density of 
material is oriented within a range of orientations near (211)f!11]. In addition, 
substantial amount of material is also oriented near (1 10)[1 12] to (1 10)1001]. One of 
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lb) 


Figure 2.8: Texture of pure Copper (a) Cu type rolling texture, (b) Cube type 
annealing texture. 
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the theories of recrystallisation says that the nucleus of recrystallisation texture 
remains in the cold rolled structure itself as small crystallites, and during the course of 
recrystallisation they just grow by consuming the highly strained region. This is the 
strain induced boundary migration theory (Block theoiy), which says that the nucleus 
of the cube texture should be present in the deformed structure itself, and it is 
observed that a small amount of material having orientations very near to the cube 
(100)[001] is also present. Figure 2.8(b) shows that most of the material has cube 
orientation for fully recrystallised alloy. The twin related component of the cube 

texture also appeared, around (221)[212]. In addition small amounts of materials are 

also found having orientation spread over from (531)[231] to (553)[5 25 ]. 

Previous experiments already proved, that among all the criteria, the presence of a 
cube oriented region in the deformed matrix is the most important one. Cube grains 
start growing from the cube oriented regions. The cube grains grow disproportionately 
larger than the grains of other orientation, because of: 

1 . Early formation, and 

2. Very high boundary migration rate. 

The grain boundary migration gives rise to bulging out of grain boundary, which is a 
very important mechanism for recrystallisation. During the recovery process, 
dislocation polygonisation gives rise to strain free subgrains, which in the later stage 
consumes the highly strained matrix in order to reduce the total free energy in the 
form of strain energy associated with dislocation and other crystal defects. 

During the experiment carried out with pure copper, a well formed subgrain having an 
orientation about 15° from the exact (001)[100] orientation has been observed to grow 
very fast consuming the neighbouring regions of other orientation (Figure 2.9). 

The extraordinary sharpness of the cube texture has been explained on the basis of 
oriented growth mechanism of recrystallisation. It is also observed that mobility of the 
grain boundary depends on the angle of the grain boundary. Low angle boundary 
shows lower mobility whereas higher angle grain boundaries exhibit higher mobility. 
In the FCC metals the fastest growing grains bear a particular orientation relationship 
to the deformed matrix. This is approximately 40° around a common <1 1 1> axis. For 
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Figure 2.9: Schematic representation of growth of Cube Texture region 

during annealing. 
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pure copper this is around 30°. Also it maintains optimum angles with Cu, Bs 
textures. Grain growth is mainly a grain boundary migration and rotation process. Due 
to the optimum angle with the neighbouring orientations, grain rotation required for 
cube orientation is low, and probably this leads to the extremely low magnitude of 
activation energy for the growth of cube region, and the cube texture dominates over 
the other orientations. 

2.10: ESTIMATION OF GRAIN BOUNDARY ENERGY DISTRIBUTION 
The properties of polycrystalline materials depend a lot on the grain boundary energy. 
The grain boundary energy state influence the strengthening properties of the 
polycrystals significantly. Several methods have been developed for the estimation of 
energy or the relative energy state of grain boundaries. Pumphery and Gleiter [Ref. 
11 ] have measured grain boundary energy by studying the spreading and 
disappearance kinetics of extrinsic grain boundary dislocations (EGBDs) in TEM 
during in-situ annealing. However, the technique is cumbersome as well as prone to 
error in the estimation of temperature and time of disappearance of EGBDs. In the 
traditional thermal groove method [Ref. 12], the original structure of the boundaries 
may undergo change. Also, the technique is limited to surface grain boundaries and 
may not be representative of the bulk material. 

Grain boundary energy may also be related to the angles between intersecting grain 
boundary planes in the bulk material. Figure 2.10 shows three grain boundaries, with 
energies yi, 72 and 73 , meeting at a triple edge. The angles (denoted as true dihedral 
angles), ai, az and as, are the angles between the grain boundary planes in the 3- 
dimensional structure of polycrystals. The true dihedral angles and the energies are 
related by the well-known Newton equation: 

= = _Il_ (1) 

sin a, sin sin aj 

Previously true dihedral angle was measured directly from the TEM images in which 
the three dimensional structure of the grain boundary planes meeting at a triple edge 
can be observed. This method of measurement is tedious and therefore it is generally 
not feasible to obtain a good statistical estimate of the distribution of true dihedral 
angles. In view of the above mentioned methods, the technique developed by Sangal 
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Figure 2 . 10 ; Sectioning of a triple edge by a plane of polish randomly oriented at 

(6, <|)). at, a2, and are the true dihedral angles and Pi, p2 and P3 are the 
corresponding plane dihedral angles; yi, 72 and 73 are the energy per unit 
area of the grain boundary planes. 
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[Ref. 13] provides an alternative approach which involves the estimation of the 
distribution of true dihedral angles from a rigorous stereological analysis of angles, 
Pi, p 2 and P 3 , between grain boundary segments (i.e., lines of intersection between the 
grain boundary planes and the sectioning plane, as shown in Figure 2.10). These 
angles are denoted as the plane dihedral angles. In a polycrystalline microstructure, 
the plane dihedral angles are the angles at the triple points, as indicated in Figure 
2 . 11 . 

Figure 2.10 shows a triple edge that is sectioned by the plane of polish P. The 
resulting linear grain boundary segments meeting at a triple point make plane dihedral 
angles. Pi, P 2 and P 3 . corresponding to the true dihedral angles, ai, 0.2 and 
respectively. It is clear from the figure that the plane dihedral angles are not only a 
function of the true dihedral angles but also depend on the orientation of the 
sectioning plane which is defined by the orientation ( 0 , ({)) of the plane normal with 
respect to the triple edge. Figure 2.12 shows that in an ideal polycrystal structure the 
triple edges are randomly oriented in 3D space and consequently the sectioning plane 
has different orientations with respect to individual triple edges. Therefore, the 
distribution of plane dihedral angles will have two components; 

( 1 ) variation in the true dihedral angles at various triple edges, and 

( 2 ) statistical variation due to the random orientation of the triple edges distributed in 
the polycrystalline structure. 

The statistical variation (like random noise) tends to overshadow or hide the actual 
distribution of true dihedral angles. Hence, such a distribution cannot be directly used 
in the estimation of grain boundary energies. However, the measurements of plane 
dihedral angles are relatively easy and it is possible to make a large number of such 
measurements in a short time. In the following section transformation procedures are 
discussed which can be used to unfold the distribution of true dihedral angles from the 
measured distribution of plane dihedral angles. 

2.11: RELATION BETWEEN PDA AND TDA; 

Figure 2.13 shows the random section of a triple edge by the plane of polish, P. The 
true dihedral angle between the grain boundary planes A and B is a and the 
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Figure 2.11: An optical microstructure of a polycrystalline material, pi, p 2 and Pa are 
the plane dihedral angles at a triple point. 
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Figure 2.12: A polycrystalline structure (X) consisting of randomly oriented triple 
edges being sectioned by the plane of polish (Ta). 
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Figure 2.13: A random section through two grain boimdary planes meeting at an 
angle of a. The normal to the plane of section makes an angle 0 with the 
z axis and the projection of the normal on the XY plane makes an angle 
with the X axis. The line of intersection of grain boundaries on plane 
of polish makes an angle p. 
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corresponding grain boundary segments a and b make a plane dihedral angle of p on 
the plane of polish as shown in Figure 2.13. The triple edge is parallel to the z-axis 
and the grain boundary plane A is in the XZ plane. The orientation of the plane of 
polish is described by the orientation of the plane normal (denoted by « ) in terms of 
angle (j) and 0 in spherical coordinates. The unit vectors and are normal to the 

grain boundary planes A and B respectively. Let a and h be the vectors parallel to 
the grain boundary segments ‘a’ and ‘b’ respectively. 

From Figure 2.13, the vectors r^, r^, and h can be expressed as follows: 

r,= J (2) 

= - sin ai + cos a j (3) 

« = sin ^ cos <j>i + sin 9 sin + cos 9k (4) 

Since a is the line of intersection of planes A and P, while h is the line of 

intersection of planes B and P, a and b can be represented by the following cross- 

products: 

a = (5) 

and 

6 = rjx« (6) 

From the equations (5) to (6), and a and b can be expressed in the following form, 
a = cos - sin ^ cos (7) 

b = cos a cos 9i + sin a cos 9j - sin 9 cos(^ - a)k (8) 

P, the plane dihedral angle, can be expressed as: 

cos p = (9) 

a.b 

From equations (7) and (8) we get: 

3*6= (cos - sin ^ cos ^i)(cos a cos 9i + sm a cos 9j -s\n9 cos(^ - a)k) 

=>3-8 = cos^ ^ cosa + sin ^ 9cos^cos(^ - a) (10) 
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a\ = -^/cos ■ ^ + sin ^ 0 cos ^ ^ = -y/l - sin ^ ^ sin ^ (j) 


( 11 ) 


= ■^1 COS' a cos^ ^ + sin ^ a cos^ ^ + sin ^ ^ cos^ - a) =.^l-sin^^sin^(^-a) ( 12 ) 


Thus from equations 2.19 to 2.22, the following expression can be obtained 

„ cos^ ^cosa + sin^ ^cos^cos(^-ar) 
cos p = , L = (13) 

y(l-sin^ ^sin^ ^)(l-sin^^sin^(^-Q:)) 

It is seen from equation 2.23 that if 0 = 0 then p = a, as would be expected. However, 
for any given value of a, P can vary from 0“ to 180“ depending upon the orientation of 
the plane of polish. 

On sectioning the polycrystalline structure (denoted as set X), the probability of an 
intersection of a triple edge (denoted as set Y) with the sectioning plane (denoted as 
set T 2 ) is given by: 

P,(T^nY\T^'t X) =—IH{X) (14) 

K 

where, / is the mean length of a triple edge, H(X) is the mean tangent diameter of X. 

Using the above, the method has been evolved for the transformation of plane 

dihedral angle distribution to true dihedral angle distribution. In this method 

individual angles at a triple edge are considered to be independent of each other. 

Consider a set Yq consisting of pairs of intersecting planes making equal angles a 

distributed in the polyctystalline space. The intersection of Ya with the sectioning 

plane, T 2 , will give rise to a probability distribution of plane dihedral angles given by 

p(P I a)dp. Let the numerical density of the set Ya be Nv(a)da, i.e., number of plane 

pairs making an angle of a per unit volume. The probability density, g(P), of plane 

dihedral angles is given by 

f i>(7’, n 7J r, t X)p{p I d)Xy ia)da 

gifi) = f T (15) 

\P(T2 nY„ \T2l-X)Ny(a)da 

<X 

Substituting P(T 2 n Ya | T 2 t X) from equation 14, the probability density, g(P) 
reduces to: 


8(P) = I P(P I a)f(a)da (16) 


29 


where, f(a) = Nv(a)y/ |^Nv(a)da is the probability density of true dihedral angles. 


Equation 16, which is an integral equation, relates the plane dihedral angle 
distribution to the true dihedral angles distribution. In order to solve equation 16, a 
numerical procedure is briefly outlined below. 

The integral equation 16, can be discretized by considering discrete distributions of 
plane as well as true dihedral angles. Let both the distributions be divided into several 
classes. The integral equation 16 is converted into a set of simultaneous equations: 


g. 

Si 


Pii P21 ••• Pii 

P12 P22 ••• Pi 2 ••• 


f. 

fl 


Sj 


Plj P2j - P.J 



(17) 


where, gj is the frequency of the class of plane dihedral angles; f, is the frequency 
of the i*'’ class of true dihedral angles; and p,j is the probability of generating a plane 
dihedral angle in the j'** class from a true dihedral angle in the i*’’ class. A Monte-Carlo 
procedure has been developed for the estimation of the coefficients py. Table 2.2 
shows the values of the coefficients p,j for various classes of a. The first and last class 
intervals have been fixed to 15° and rests all were 10°. For different class widths, a 
new set of values for these coefficients has to be calculated. The frequency gj can be 
determined by the measurement of plane dihedral angles on the polycrystalline 
microstructures. Now the frequency distribution of true dihedral angles can be 
calculated by solving the simultaneous equations given in relation 17 [Ref. 14]. 
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Table 2.2. Computed values of the coefficients a, 


PDA distribution (a,j) for different value of TDA 
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CHAPTERS 

MATERIAL, MICROSCOPY AND IMAGE ANALYSIS 


3.1: MATERIAL AND INITIAL TREATMENT 


The chemical composition (both nominal as well as the detailed compositions) of the 
alloys (wt%) used in the present study are shown in Table 3.1. 

Table 3.1: Chemical Composition of the Alloys 


Alloy Nominal 

Composition 

Detailed Weight Percent of the Alloying Elements 

C 

Si 

Mn 


S 

Ni 

Co 

Ni-10%Co 

0.006 

<0.03 

<0.03 

0.005 

0.003 

88.40 

11.15 

Ni-20% Co 

0.007 

0.060 

<0.03 

0.005 

0.003 

76.90 

22.85 

Ni-30% Co 

0.006 

<0.03 

<0.03 

0.004 

0.003 

68.85 

30.90 

Ni-40% Co 

0.006 

0.030 

<0.03 

0.003 

0.004 

58.70 

41.05 

Ni-60%Co 

0.006 

0.060 

<0.03 

0.003 

0.004 

1 

39.20 

60.50 


All the alloys were induction melted and cast under an argon atmosphere using virgin 
metals of purity greater than 99.99%. The ingots were cylindrical in shape having a 
height of about 20 mm. These ingots were cold rolled 50% to a thickness of 10 mm 
and then homogenisation annealed in vacuum at 1150°C for 24 hours. These were 
then again cold rolled 50% and annealed at 1 150°C for a period of three hours to yield 
the starting material of almost random texture and a grain size of ~ 0.1 mm. 

3.2: COLD ROLLING 

Strips cut from the thick sheets were cold rolled in 2-high laboratory rolling mill 
where paraffin oil was used as lubricant. The ratio of the length of contact with the 
rolls to specimen thickness was maintained at a value greater than unity in order to get 
homogeneous deformation throughout the specimens. During rolling the strips were 
reversed after each pass and between two successive passes the strips were immersed 
in cold water to reduce the temperature rise due to deformation. Each alloy was rolled 
to 95% final reduction in thickness. 

3.3: ANNEALING 

After cold rolling the samples were annealed. Annealing was carried out in a 
Laboratory Tube type furnace. Each of the five samples of square shape was 
suspended with the help of wires inside the furnace. All samples were annealed at 












































800°C for seven different timings of 15 mins, 30 mins, 1 hour, 3 hour, 10 hours, 20 
hours and 50 hours respectively in an Argon atmosphere. Gas was continuously 
passed through the furnace, and finally the samples were quenched in normal water. 
3.4: TEXTURE MEASUREMENT 

Texture measurements were made on specimens taken from all the recrystallised 
samples. Specimens of the si/e (25 mm x 1 5 mm) were cut from the recrystallised 
material for this purpose. All the texture specimens were lacquered on one face and 
half of the thickness was removed from each of them chemically in a solution of 1 1 % 
HF and 89% HNO 3 . All measurements were done on the mid-thickness plane of each 
sample. Four pole figures {111}, {200}, {220} and {311} were measured by the 
Schultz reflection method [Ref. 15] for each sample. From these data the ODF 
(Orientation distribution function) was determined, using the method of Bunge [Ref. 
9]. Finally, modelling of each ODF was carried out using Gauss type functions for the 
different components and the volume fraction of each component was determined 
there from. [Ref. 16]. 

The texture measurements were carried out in the Institut fur Metallkunde, Aachen in 
Germany. 

3.5: OPTICAL MICROSCOPY 

The cold rolled samples were cold mounted for examination under the optical 
microscope. The samples for optical microscopy were ground and polished on emery 
papers (0 to 4 grades) followed by wet polish with alumina suspension of particle size 
1 pm and 0.03pm respectively. These were subsequently etched before examining the 
microstructure. The etching reagents used are: 

(a) 5 gm Ammonium Persulphate, 25 ml H 2 O. 

(b) 83 gm FeCls, 33 ml. HCl. 

(c) lOml.HNOs. 

Microstructures were observed on Lieca Optical Microscope for all the samples. 

3.6: IMAGE PROCESSING ALGORITHM 

24 bit RGB digital images were acquired through a CCD camera mounted on Lieca 
Optical Microscope. These images were subsequently converted to 8 -bit grayscale 
images. The size of the acquired image was 640X480 picture elements. A micro-scale 
was used to perform calibration of the image. At a magnification of lOOX, dimension 


33 



of one pixel was calculated to be 0.6041 jm . The software used was Image-Pro Plus 
ver 4. 1 for capturing as well as processing the images. 

3.6.1; BINARIZATION/ THRESHOLDING 

A typical grayscale image of the microstructure of Ni-20%Co and the corresponding 
bimodal gray level distribution are shown in Figure 3.1(a) and 3.1(b) respectively. 
For the binarization of such images, automatic thresholding was applied, which 
involved fixing the threshold limit at a minimum between two maxima of gray level 
distribution. All pixel values below the threshold limit were set to ‘0’ value and those 
pixels whose value was above this limit were fixed to ‘255’ value. This type of 
binarisation gave consistent results independent of image contrast and brightness 
[Ref. 17]. Figure 3.2 shows the resulting binary image in which all the sensor noise 
and background illumination variation have been abstracted. 

3.6.2: DILATION FILTERING 

The dilation filter is a morphological filter that changes the size of an object in an 
image by dilating the specified objects. Actually it dilates bright objects, and erodes 
dark ones. The areas, which had to clean up, were selected as shown in Figure 3.2 
and then dilation filter was applied. Some narrow peninsulas or small particles are 
disappeared after dilation. This filter remove noise/artifacts created by etch pits, 
inhomogeneous etching/polishing etc. in the image as shown in Figure 3.3. 

3.6.3: EROSION FILTERING 

The erosion filter is also morphological filter that changes the size and shape of an 
object in an image by eroding or reducing the boundaries of the specified objects. It is 
a shrinking filter therefore eroded objects are smaller than the initial one. Basically 
the erosion filter erodes the edges of bright objects, and enlarges dark ones. After 
using this filter the discontinuity of the grain boundary is removed. This operation 
also resulted in the removal of some narrow peninsulas and small artifacts as shown in 
Figure 3.4. 

3.6.4: THINNING MORPHOLOGICAL FILTER 

The eroded image is transformed by a thinning morphological operation to reduce the 
grain boundaries to only one pixel thick line. The thinning filter was applied 
iteratively till further operations of the filter had no effect on the image. The resulting 
thinned image is shown in Figure 3.5. 
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Figure 3.1: (a) A typical grayscale image of a polycrystalline microstracture, 
(b) Bimodal gray level histogram of the image (a). 
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Figure 3.2: Transformed binary image (after applying thresholding algorithm) 
showing polishing/etching artifacts. 
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Figure 3.3: Result of applying dilation filter on Figure 3.2. The disappearance of 
peninsulas and etching artifacts is noticeable. 
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Figure 3.4; Image of microstructure obtained after erosion filter was applied on the 
image shown in Figure 3.3. It is noticeable here that the grain boundary has become 

continuous. 
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Figure 3.5; Image of the resultant microstructure after application of thinning filter 
on the image shown in Figure 3.4. The grain boimdaries are now only one pixel thick. 
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3.6.5; BITMAPS 

The final image was converted into a bitmap. Figure 3.6 shows a bitmap of the 
rectangular region marked in Figure 3.5. As can be seen from this figure, the grain 
boundary pixels are represented by O’s and the matrix pixels are represented by 1 ’s 
(these pixels uniquely had a value of 255 in the binary image). The bitmap was used 
as input to the algorithm developed to detect triple junctions and determine their 
geometry, as discussed below. 

3.7: TRIPLE POINT DETECTION 

Refer to the bitmap of the rectangular section of the image (Figure 3.5) of 
microstructure as seen in Figure 3.6. Each grain boundary pixel (pixel value of ‘0’) 
was considered at the center of a 3x3 matrix. If a grain boundary pixel had exactly 
three other boundary pixels in its neighborhood then the pixel was marked as a 
possible triple point, and subjected to further evaluation by comparing it with a pre- 
generated list of valid triple point configurations. Partial lists of valid configurations 
in the form of 3x3 matrices are shown in Figure 3.7. The remaining configurations 
can be generated by rotating the individual matrices through 3t/2, n and 3it/2 thereby 
generating a total of 12 valid configurations. However there are many other 
configurations possible in which the central ‘0’ pixel has three neighboring pixels, but 
are not included in the above list. The reason for this becomes evident from 
Figure 3.8. The (row, column) coordinates of each pixel in the bitmap can be referred 
relative to the left-top pixel, whose coordinate is assigned (1,1), as shown in Figure 
3.8. On examining this figure, it is clear that grain boundaiy pixels at coordinates 
(3,6), (3,7), (4,7) and (4,8) are not triple points even though each of these pixels have 
exactly three neighboring grain boundary pixels, such triple points have been termed 
as pseudo triple points [Ref. 18]. Similarly the grain boundary pixels at coordinates 
(3,4) and (4,3) are also pseudo triple points. However the triple point at coordinate 
(3,3) is an actual triple point, which matches with the configuration shown in Figure 
3.7a rotated through /r . Thus, the pseudo triple points can be classified according to 
the following general rule: If any two of the three neighboring pixels share a boundaiy 
among themselves, which means that if any two neighboring pixels have same ‘x’ 
(ordinate and the ‘y’ coordinate differs by one and vice-versa then the triple point is 
assified as a pseudo triple point. This rule was found to be appropriate by comparing 
e results of the automated detection of triple points with manual detection. 
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Figure 3.6: Bitmap of the marked rectangular region in Figure 3.5. ‘0’ refers to grain 
boundary pixels and ‘1’ refers to matrix pixels. 
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Figure 3.7: A partial configurations of valid triple points. Remaining configurations 


can be generated rotating the matrices through 7c/2, % and 37c/2. 
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Figure 3.8: Illustration of pseudo triple points. 
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3.8: EVALUATION OF PLANE DIHEDRAL ANGLES (PDA) 

After detecting triple points, as described in section 3.7, the coordinates of the pixels 
of the grain boundary segments were traced till the segment terminated into another 
triple point or till no further touching pixels could be found. The algorithm developed 
for the same is at appendix A.l. At this stage, if a particular segment was found to be 
very short (less than four pixel length) then the associated triple point was considered 
as a pseudo-triple point and therefore eliminated from the analysis. Manual 
examination of the microstructure confirmed the validity of this procedure. 

A least square procedure was employed to fit polynomials on the pixels of the grain 
boundary segments at each triple point. The form of the polynomial chosen for the 
purpose is given as: 

y = ao+a,x + a 2 X^ +... + a„x'’' (3.1) 

where ‘m’ is the polynomial order, ai, a 2 am are the coefficients. This form of the 
polynomial was modified to ensure that the fitted curve is constrained to pass through 
the triple point. 

A least square procedure was adopted to determine the values of the polynomial 
coefficients. This involved minimization of the sum of square (R^) of the following 
equation: 




/ = l 


(3.2) 


where, y, is the measured value and y is given by equation (3.1) for the ith pixel. 
Derivative of equation (3.2) with respect to each of the unknown coefficients of the 
polynomial is performed and equated to zero to develop a set of equations, which was 
solved by gauss-elimination technique. 

Each grain boundary segment was fitted with polynomial of increasing degree of 
order, starting with the order of 1 . In order to select the best polynomial order, a ratio 
F was evaluated to compare the goodness-of-fit (quantified by sum of squares, R^) of 
two polynomials of orders m and m+1 



(3Ja) 


and 
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y^=n-{m + \) and y 2 =n-{m + 2) (3.3b) 

where, degree of freedom, n = number of data points and the subscripts 1 and 2 

refer to the polynomials of orders m and m+1 respectively. 

If the fit of the higher order polynomial (order = m+1) has a larger sum of squares, 
i.e., R 2 ^>Ri^, then clearly the fit becomes worse on increasing the order and therefore, 
the lower order polynomial (order = m) is chosen as the best order. In terms of the F- 
ratio (equation 3.3a), the value of F will be negative. However, in general with 
increasing polynomial order the sum of square value is expected to reduce because the 
higher order equation has more adjustable parameters. Thus, a further analysis of the 
F-ratio was done. 


In going from lower order to the higher order polynomial, if the relative decrease in 
the sum-of-squares was less than the relative decrease in the degree of freedom (i.e., 
F<=1), then the polynomial with the lower order was selected as the best polynomial 
order. On the other hand, if the value of F was greater than I, i.e., the relative 
decrease in the sum of squares was larger than the relative decrease in the degree of 
freedom, then the higher polynomial order was accepted and the above analysis was 
repeated by further incrementing the polynomial order. 

The correlation coefficient, r, for each polynomial fit was determined by using the 
following equation: 



where, . 

/*=! 

For a perfect fit R^ = 0 and r^ = 1, signifying that the polynomial fit truly fits all the 
coordinates of each grain boundary. While close to zero signifies a poor fit. The 
correlation coefficient was used to decide whether to accept or reject a particular fit, 
as discussed in the section. 


The plane dihedral angles were determined as the angles between the tangents to the 
grain boundary segments at the triple points. Figure 3.9 illustrates a typical example 
of three tangent lines (calculated by differentiating the fitted curves) at a triple point. 
Plsme dihedral angle, P, is given by 


cos>0 = 


(pAf +{OBf -{ABf 
2iOA)(OB) 


(3.5) 
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Figure 3.9: Tangents (thick lines) drawn to the fitted grain boundary segments 

(thin lines) at a triple point ‘O’. ‘A’ and ‘B’ are arbitrary points on the 
segments. 




Figure 3.10: Schematic of cumulative distribution for true dihedral angles: 

(a) Omax » (b)atnm. Pi and p 2 are randomly generated number, xj 
and \2 are corresponding Omax and amm 



Similarly, the other two plane dihedral angles can also be determined. 

After calculating the plane dihedral angles of each triple point the maximum and 
minimum plane dihedral angles of each triple points were evaluated. To find out the 
distribution of plane dihedral angles 17 small classes were considered over an angular 
range of 0 ° to 180 °. The first and last class intervals were 15 ° and rests all were 10 °. 
Then the frequency of each class was calculated for both maximum and minimum 
plane dihedral angles. Subsequently the relative frequencies were computed for both. 
3 . 9 : TRANSFORMATION OF PDA DISTRIBUTION TO TDA DISTRIBUTION 
The relative frequencies of plane dihedral angles (both maximum and minimum) were 
transformed to true dihedral angles by the methods as described in section 2.1 1 with 
the help of computer program. The minimum and maximum relative fi-equencies of 
true dihedral angles were denoted as Omm and Omax respectively. For transformation the 
values of coefficients a,j were used as mentioned in Table 2 . 2 . 

3 . 10 : RELATIVE GRAIN BOUNDARY ENERGIES 

The relative grain boundary energy was calculated in terms of two parameters: 

= ( 3 . 6 ) 

^3 sin 0:3 


^2 _sina2 

^3 sin 


( 3 . 7 ) 


where yi< 72 < 73 are the energies of grain boundary planes meeting at a triple edge. In 
order to calculate above, joint distribution of Omm, Omax was computed using a Monte- 


Carlo procedure. For this the calculated cumulative distributions of amin» Omax were 
taken, as schematically shown in Figure 3 . 10 (a) and (b). 

The Monte-Carlo procedure involved generation of a pair of uniformly distributed 
random numbers (in the range: 0 to 1 ) at each iteration. Let a pair of such numbers be 
pi and p2. For each pi and p2 the corresponding xi and X2 could easily be found out. 
71, 72 and 73 were then computed from equation 3.6 and 3 . 7 . 5000 iteration was done 
fi-om which we obtained the distribution of both 71/73 and 72/73. The relative 
frequencies of 71/73 and 72/73 were evaluated subsequently. The whole process was 
done with the help of computer program as mention in appendix A. 7 . 



CHAPTER 4 

RESULTS AND ANALYSIS 


As mentioned in chapter 3, the five Ni-Co alloys i.e. Ni-10%Co, Ni-20%Co, Ni- 
30%Co, Ni-40% Co and Ni-60%Co were cold rolled by an amount of 95% reduction 
to approximately 0.25 mm thickness. In ord'er that homogeneous deformation through 
out the specimens might be achieved, the ratio of the length of contact with rolls to 
specimen thickness was maintained at a value greater than unity. To produce 
rccrystallisation texture the cold rolled samples were annealed at 800”C for seven 
different times of 15 minutes, 30 minutes, 1 hour, 3 hours, 10 hours, 20 hours and 50 
hours respectively. 

4.1: OPTICAL MICROSCOPY 

The optical microstructures of annealed samples of all the five alloys are shown in 
Figure 4.1 to Figure 4.7. In general all the samples showed almost full 
recrystallisation from the shortest annealing time of 15 minutes. Grain growth was 
also observed in the samples. Annealing twins can also be observed in almost all the 
samples. Each microstructure as developed in the normal manner was further 
processed with the help of an Image Analyser as discussed in chapter 3. Some typical 
processed microstructures are shown in Figure 4.8 and 4.9. The processed 
microstructures appear far more clear as compared to the unprocessed microstructures 
indicating very clearly the grain boundaries and triple points. The important thing to 
be noticed is that the twin boundary areas were eliminated from the processed 
structure. The polishing/etching artefacts were also removed. After processing the 
grain boundary regions were one pixel thick which made it easier to determine the 
plane dihedral angles. For each sample nearly 150 triple points were considered from 
which plane dihedral angles were determined. The plane dihedral angle distributions 
were converted to true dihedral angle distributions and subsequently relative grain 
boundary energy ratios were evaluated, in the manner described in chapter3. 

4.2: RELATIVE GRAIN BOUNDARY ENERGY 

The relative grain boundary energy ratios, both yl/y3 and y2y^3 were calculated by the 
method described in chapter 3 where yl<y2<y3. The frequency distributions of these 
values have been plotted against the values of these ratios for all the samples 
concerned, in different manners to arrive at meaningful conclusions. 
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(c) Ni-30% Co 


(d)Ni-40%Co 



(e)Ni-60%Co 


Figure 4.1: Micrographs of different Ni-Co alloys annealed for 15 minutes. 
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(a)Ni-10%Co 


(b)Ni-20%Co 



(c) Ni-30% Co 


(d)Ni-40%Co 



(e)Ni-60% Co 


Figure 4.2: Micrographs of different Ni-Co alloys annealed for 30 minutes. 









(c)Ni-30%Co 


(d)Ni-40%Co 



Figure 4.4: Micrographs of different Ni-Co alloys annealed for 3 hours. 








(c) Ni-30% Co 


(d)Ni-40%Co 



(e)Ni-60%Co 

Figure 4.5: Micrographs of different Ni-Co alloys annealed for 10 hours. 








(e)Ni-60%Co 


Figure 4.6; Micrographs of different Ni-Co alloys annealed for 20 hours. 
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(b)Ni-20%Co 



(c)Ni-30%Co 


(d)Ni-40%Co 



Figure 4.7; Micrographs of different Ni-Co alloys annealed for 50 hours. 








(c)Ni-30%Co (d)Ni-40%Co 



(e)Ni-60%Co 


Figure 4.8: Micrographs showing processed structure of different Ni-Co alloys annealed 

for 1 hour. 





(c) Ni-30% Co 


(d)Ni-40%Co 



(e)Ni-60%Co 


Figure 4.9: Micrographs showing processed structure of different Ni-Co alloy annealed 

for 50 hours. 








4.2.1: HISTOGRAMS 


4.2.1.1: HISTOGRAMS OF FREQUENCY vs. yl/yS 

The histograms of frequency vs. relative grain boundary energy ratio yj/ys are shown 
in Figures 4.10 to 4.16. It is observed from the histograms that the relative grain 
boundary energy ratio with a class width of 0.6 to 0.7 and 0.7 to 0.8 always possess 
higher frequency compared to others. Relative grain boundary energy ratios with a 
value greater than 0.9 are never obtained in any of the samples. 

4.2.1.2: HISTOGRAMS OF FREQUENCY vs. y2/y3 

'fhe histograms of frequency vs. relative grain boundary energy ratio yihi are shown 
in Figure 4.17 to 4.23. In contrast to the histograms of frequency vs. yi/ys, in this case 
the y 2 /y 3 ratios always possess values higher than 0.6. The frequency of ya/ys ratios 
with a value less than 0.6 is practically zero. This has been observed in all the 
samples. To be precise, the ya/ys ratios have minimum frequencies in the class width 
of 0.6 to 0.7. That means ya/ys ratios of values higher than 0.7 have higher frequencies 
compared to others. The frequency of ya/ys ratio with a class \vidth of 0.9 to 1 is 
maximum for majority of the cases. 

4.2.2; FREQUENCY DISTRIBUTION OF RELATIVE GRAIN BOUNDARY 
ENERGY RATIO 

All the relevant data that were used to determine the histograms were plotted in a 
graphical form for simplicity and for ease in comparison. 

4.2.2.1: FREQUENCY DISTRIBUTION OF yl/y3 

The frequency distribution of yi/yj ratio for various annealing times are shown in 
Figure 4.24 to 4.28. Upto certain value of yi/ys ratio, there is no significant difference 
in frequency. This trend changes from around 0.2 (value of yi/ys ratio), after which 
there is a continuous increase in frequency. The peaks are distributed in the region 0.6 
to 0.8. After attaining the peak value there is a drastic decrease in frequency and at 
yi/y 3 ratio of 1 the frequency is zero. This is quite obvious because yi and yi will 
never have the same value except under perfect equilibrium condition. 

4.2.2.2: FREQUENCY DISTRIBUTION OFy2/y3 

The frequency distribution of ya/ya ratio for various annealing time are shown in 
Figure 4.29 to 4.33. For all the alloys the graphs show similar trends. Upto the ya/ys 
ratio of 0.6 the frequency is zero irrespective of different alloys and annealing times. 
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Figure 4.11; Frequency distribution of relative energy ratio yl/yS of different Ni-Co 

alloys annealed for 30 minutes. 
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(e)Ni-60%Co 


Figure 4.13; Frequency distribution of relative energy ratio yl/y3 of different Ni-Co 

alloys annealed for 3 hours. 
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Figure 4.14: 


Frequency i Frequency 


22-1 



(b)Ni-20%Co 


24-1 



rl/y3 

(d)Ni-40%Co 



0 6 0.7 0.S 0.9 1 










Frequency frequency 








Frequency 



(c) Ni-30% Co 


(d)Ni-40%Co 


05 
0.45 
04 
0.35 
c 03 
3 025 

f 02 

it 015 
01 
005 
0 


01 020304050.6070600 1 


(e)Ni-60%Co 

Figure 4 . 16 : Frequency distribution of relative energy ratio yl/yS of different Ni-Co 

alloys annealed for 50 hours. 
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Figure 4.17; Frequency distribution of relative energy ratio yUyZ of different Ni-Co 

alloys annealed for 15 minutes. 
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(c) Ni-30% Co (d) Ni-40% Co 



(e) Ni-60% Co 

Figure 18: Frequency distribution of relative energy ratio Y2/y3 of different Ni-Co 

alloys annealed for 30 minutes. 
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(c) Ni-30% Co (d) Ni-40% Co 
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Figure 4.21; Frequency distribution of relative energy ratio y2/y3 of different Ni-Co 

alloys annealed for 10 hours. 
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Figure 4.26: Frequency distribution of relative grain boundary energy ratio yl/y3 for 
Ni-30% Co samples subjected to different annealing times. 
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Figure 4.27; Frequency distribution of relative grain boundary energy ratio yl/yS for 
Ni-40% Co samples subjected to different annealing times. 
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Figure 4.28: Frequency distribution of relative grain boundary energy ratio yl/y3 for 
Ni-60% Co samples subjected to different annealing times. 







Figure 4.29: Frequency distribution of relative grain boundary energy ratio iHyZ for 
Ni-10% Co sanqjles subjected to different annealing times. 
















Figure 4.32; Frequency distribution of relative grain boundary energy ratio Y2/y3 for 
Ni-40% Co samples subjected to different annealing times. 
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Figure 4.33: Frequency distribution of relative grain boundary energy ratio 72/73 for 
Ni-60% Co samples subjected to different annealing times. 





Beyond that there is an increase in frequency upto the 72/73 ratio of 1. The fact is that 
the distribution of grain boundary energies is quite similar and this behavior is 
maintained for all the alloys. Although there is some small variation in the plots, it is 
not significant. 

4.2.3: VARIATION OF FREQUENCIES WITH ANNEALING TIME 
4.2.3.1: FOR RELATIVE ENERGY RATIO 7 I /73 

The frequencies of different classes of 7 I /73 as a function of annealing time in the 
five Ni-Co alloys are shown in Figures 4.34 to 4.38. 

For Ni-10% Co alloy after 20 hours of annealing the frequencies of different classes 
ofyl/yS become stabilized. 

For Ni-20% Co the classes 0.5 to 0.6 and 0.7 to 0.8 of 7 I /73 increases drastically after 
20 hours of annealing at the expense of other components. 

for Ni-30% Co after 20 hours of annealing the 0.6 to 0.7 class of 7 I /73 starts growing 
drastically at the cost of other components. 

For Ni-40% Co and Ni-60% Co all other classes of 7 I /73 except 0.6 to 0.7 and 0.7 to 

0. 8 reach an equilibrium value after 20 hours of annealing. From the initial stage of 
annealing they grow at the expense of each other. After 50 hours of annealing they are 
present in equal frequencies. 

4.2.3.2: FOR RELATIVE ENERGY RATIO 72/73 

The frequencies of different classes of 72/73 as a function of annealing time in the 
four Ni-I'e alloys are shown in Figures 4.39 to 4.43. In this case the important classes 
of 72/73 are 0.7 to 0.8, 0.8 to 0.9 and 0.9 to 1. 

In Ni-10% Co alloy the classes 0.8 to 0.9 and 0.9 to 1 of 72/73 grow at the expense of 
one another. 

The class of 72/73 having maximum frequency in Ni-20% Co and Ni-30% Co is 0.9 to 

1 . 

In Ni-40% Co and Ni-60% Co alloys the classes 0.8 to 0.9 and 0.9 to I of 72/73 grow 
at the expense of one another. 

4.2.4: AVERAGE vs. ANNEALING TIME 

The variations of average values of 7l/73and 72/73 for the different alloys as a 
function of annealing time are shown in Figures 4.44 and 4.45 respectively. After 20 
hours of annealing Ni-10% Co alloy showed highest average value for y^h(^ which 
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Figure 434; The variation of frequencies of different classes of yl/y3 as a function of 

annealing time for Ni-10% Co samples. 
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Figure 435: The variation of frequencies of different classes of yl/y3 as a frinction of 

annealing time for Ni-20% Co samples. 
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Figure 4.36: The variation of frequencies of different classes of yl/yS as a function of 

annealing time for Ni-30% Co samples. 
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Figure 4.37: The variation of frequencies of different classes of yl/yS as a function of 

annealing time for Ni-40% Co samples. 
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Figure 4.38: The variation of frequencies of different classes of yl/yS as a function of 

annealing time for Ni-60% Co samples. 
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Figure 4.39: The variation of frequencies of different classes ofy2/y3 as a function of 

annealing time for Ni-10% Co samples. 
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Figure 4.40: The variation of frequencies of different classes of y2/y3 as a fiinction of 

annealing time for Ni-20% Co samples. 
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Figure 4.41: The variation of frequencies of different classes of y2/y3 as a function of 

annealing time for Ni-30% Co samples. 
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Figure 4.42: The variation of frequencies of different classes of y2/y3 as a function of 

annealing time for Ni-40% Co samples. 
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Figure 4.43; The variation of frequencies of different classes ofyUyZ as a function of 

annealing time for Ni-60% Co samples. 
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Figure 4.43: The variation of frequencies of different classes of yl/yS as a function of 

annealing time for Ni-60% Co samples. 
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Figure 4.44: The variation of mean of relative energy ratio yl/yS of different alloys as 

a function of annealing time. 
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Figure 4 . 45 : The variation of mean of relative energy ratio 72/73 of different alloys as 

a function of annealing time. 



indicates that it has reached more near towards equilibrium than the other alloys. In 
the r2/Y3 plots there is hardly any variation among the alloys. 

4.2.5: COEFFICIENT OF VARIATION vs. ANNEALING TIME 
For the purpo.se of comparison, the coefficient of variation (i.e., standard 
deviation/mean) of yl/yS and y2/y3 as a function of annealing time are plotted in 
Figure 4.46 and 4.47 respectively. In case of yl/y3 there is a significant scatter of 
coefficient of variation. In contrast to yl/y3, the fluctuation of coefficient of variation 
is less for y2/y3. This is quite expected, as the dihedral angles become closer to each 
other and the grain boundary energy variation also reduces. 

4.3; TEXTURE RESULTS 

The various texture components in different alloys for different aimealing times are 
given in Tables 4.1 to 4.30. From the tables it can be observed that in the samples of 
Ni-10% Co, Ni-20% Co, and Ni-30% Co the cube component {100}<001> is 
predominantly present. Along with cube component its twin components are also 
present. The first generation twin of cube is represented as CTl {212}<122> and 
second generation twins of cube are represented as CT2 {148}<841> and 
{148}<474>. In Ni-40% Co alloy, cube component is present in little amount. In Ni- 
10% Co, Ni-20% Co, and Ni-30% Co alloys a component {437}<4 11 7> which is 
near to first generation twin of Bs recrystallisation component (BsT) is also present 
consistently. In Ni-40% Co alloy three components are important: Bs recrystallisation 
component (BsR) (236}<385>, {279}<7 11 7>, and the retained S component 
{123}<634>. In Ni-60% Co alloy the primary components are rotated Goss 
component {01 1 }<1 10> and a component which is near to it i.e. {01 1 }<1 14>. 

4.3.1: VARIATION OF CUBE COMPONENT WITH ANNEALING TIME IN 
Ni-10% Co TO Ni-40% Co ALLOY 

The variations of intensity and volume fiaction of cube component as a fiinction of 
annealing time for different alloys are shown in Figures 4.48 and 4.49 respectively. 
The intensity of cube component in both Ni-30% Co and Ni-40% Co alloys do not 
vary much and remains nearly the same for almost all annealing times. Although Ni- 
10% Co and Ni-20% Co show larger variation at the initial period of annealing time, 
after 20 hours of annealing the intensity of all the allo)^ become steady. In case of 
volume fraction of cube conqwnent upto 10 hours of annealing time there are some 
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Figure 4.46: The variation of co-efficient of variation of relative energy ratio yl/y3 of 
different alloys as a function of annealing time. 
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Figure 4.47: The variation of co-efficient of variation of relative energy ratio y2/Y3 of 
different alloys as a function of annealing time. 




Table: 4.1: Ni-10% Co annealed at 800°C for 15 minutes. 


Serial No. 


Component 

ID 


Rotated 

Cube 


Cube 


CT1 




0 


66.6 


Suler Angles 


Jl 


34.4 


0.4 


45.8 


(b 


62.6 


f(g) 


3.53 


19.25 


6.76 


Volume % 


6.59 


29.36 


30.08 


Near to Bs 
Rex. T1 


CT2 


61.2 

74.4 


35.1 


50 


28.9 


75 


2.6 


2.3 


14.7 


10.55 


6 


CT2 


14.7 


25.7 


10 


2.3 


7.95 


Table: 4.2: Ni-10% Co annealed at 800°C for 1 hr. 



Component 

ID 


Euler Angles 


1 

Rotated 

Cube 

0 

30.7 

2 

Cube 

0 

0.4 

3 

CTl 

65 

44.6 

4 

Near to Bs 
Rex. T1 

61.2 

36 

5 

CT2 

15.3 

25.2 

6 

CT2 

71.1 

29.1 


h 

f(g) 

Volume % 

0 

4.1 


0 

23.22 

41.62 

62.8 

5.62 

21.22 

55 

5.1 

18 

10 

2.4 

5.33 

75 

2 

5.76 


Table: 4.3: Ni-10% Co annealed at 800°C for 3 hr. 


Serial No. 

Component 

ID 

] 

Euler Angles 

f(g) 

Volume % 


|||||||||!||[||||^^ 

h 

1 

Rotated 

Cube 

0.5 

25.9 

0 

4.36 

5.11 

2 

Cube 

0 

0.4 

0 

32.59 

45.66 

3 

CTl 

64.2 

45.2 

61.8 

7.02 

18.53 

4 

Near to Bs 
Rex. T1 

60.8 

i 

35.8 

55 

6 

14.41 

5 

CT2 

71.4 1 

30.6 

75 

2.3 

9.78 

6 

CT2 

15.3 

■ 25.5 

10 

2.6 

3.59 








































































































Table: 4.4: Ni-10% Co annealed at SOOT for 10 hr. 



Component 

ID 

1 

Euler Angles 


Volume % 




1 

Rotated 

Cube 

0 

25.4 

0 

5.12 

5.09 

■ 2 

Cube 

0 

0.5 

0 

27.76 

42.81 

3 

Near Bs 
Rex. T1 

58.3 

35.8 

57.5 

4.98 

21.46 

4 

CTl 

65.3 

44.4 

60 

5.33 

20.4 

5 

CT2 

16.4 

25.7 

10 

2.5 

3.79 

6 

CT2 

74.2 

30.1 

75 

2.2 

3.93 


Table: 4.5: Ni-10% Co annealed at 800°C for 20 hr. 


Serial No. 

Component 

ID 

1 

Euler Angles 

f(g) 

Volume % 




1 

Cube 

0 

0.4 

0 

34.56 

44.71 

2 

CT2 

15.7 

26.7 

11.9 

2.88 

9.96 

3 

Near to Bs 
Rex. T1 

60.4 

35.1 

55.6 

5.84 

22.57 

4 

CTl 

65 

45 

61.8 

7.07 

14.1 

5 

Goss 

0 

45 

0 

4.5 

1.48 

6 

CT2 

70.8 

29.9 

75 

2.6 

4 

7 

Rotated 

Cube 

0 

28.5 

0 

4.5 

2.06 


Table: 4.6: Ni-10% Co annealed at 800°C for 50 hr. 


Serial No. 


Component 

ID 


Rotated 

Cube 


Cube 


Goss 


CTl 


CT2 


Near to Bs 
Rex. T1 


CT2 


Euler Angles 

f(g) 

Volume % 





25.1 

0 

5.8 

6.06 

0.4 

0 

32.73 

46.91 

44 

0 

4.69 

3.24 

44.6 

61.2 

6.67 

14.2 . 

25.3 

10 

2.5 

3.58 

35.1 

50 


15.94 


71.2 


29.5 


75 


2.5 


7.94 



























































































































Table: 4.7: Ni-20% Co annealed at 800°C for 15 minutes. 




















































































































Table: 4.10: Ni-20% Co annealed at 800®C for 10 hr. 



Table: 4.11 : Ni-20% Co annealed at 800°C for 20 hr. 



















































































































Table: 4.13: Ni-30% Co annealed at 800T for 15 minutes. 


Serial No. 

Component 

ID 

] 

Euler Angles 

f(g) 

Volume % 


(t> 

<l>2 

1 

Rotated 

Cube 

0 

30.4 

0 

3.13 

7.4 

2 ' 

CT2 

16 

25.9 

3.6 

2.28 

15.38 

3 

Cube 

0 

0.4 

0 

9.53 

15.94 

4 

{019}<991> 

45 

0 

0 

0.56 

3.54 

5 

Near to Bs 
Rex. T1 

61.1 

35.6 

55 

3.26 

15.81 

6 

CTl 

65.8 


63.2 

3.9 

10.73 

7 

CT2 

71.2 

31 

70 

2.34 

22.75 

8 

{416}<121> 

44.8 

39.8 

75 

1.94 

r 5.86 


Table: 4.14: Ni-30% Co annealed at 800°C for 1 hr. 


Serial No. 



Component 

ID 


{0 611}<0 
11 6 > 


Rotated 

Cube 


Goss 


Cube 


CT2 


Near to Bs 
Rex. T1 


Bs Rex. T1 


CT2 


CTl 


Euler Angles 


Volume % 



10.58 


.2 


3.99 








































































































































Table: 4.16: Ni-30% Co annealed at SOO^C for 10 hr. 



Component 

ID 

] 

Euler Angles 

f(g) 

Volume % 

miiinsin 


<1>2 

1 

Rotated 

Cube 

0 

— _ — 

26.3 

0 

3.17 

2.95 

2 

Cube 

■HHI 

0.5 

0 

15.31 

23.4 

3 

Goss 

HKHI 

45 

0 

2.35 

1.79 

4 

Near to Bs 
Rex. T1 

59 

35 

57.3 

6.69 

42.39 

5 

CT2 

75.4 

30.5 

70 

2 

14.4 

6 

CT2 

14.6 

25.9 

10 

1.6 

6.6 

7 

CTl 

55 

40 

65 

5.2 

7.56 


Table: 4.17: Ni-30% Co annealed at 800®C for 20 hr. 


Serial No. 

Component 

ID 

] 

Euler Angles 

%) 

Volume % 


llllllllllllllllll^^ 

<l>2 

1 

Rotated 

Cube 

0 

25.3 

0 

3.2 

3.99 

2 

Goss 

0 

45 

0 

3.31 

1.66 

3 

Cube 

0 

0.5 

0 

19.84 

25.72 

4 

Near to Bs 
Rex. T1 

61.6 

35 

54.6 

8.01 

33.51 

5 

CTl 

60 

45 

65 

4.7 

5.57 

6 

CT2 

70 

30 

75 

2.1 

23.2 

7 

CT2 

16 

26.7 

10 

1.9 

4.85 


Table: 4.18: Ni-30% Co annealed at 800“C for 50 hr. 


Serial No. 



Component 

ID 


Rotated 

Cube 


Cube 


Goss 


Near to Bs 
Rex. T1 


CTl 


CT2 


CT2 


Euler Angles 


Volume % 



60 

14,5 

74.3 


45 

25.6 

30 


65 

jO 

75 


M 

2.6 


10.53 

3.23 

8.31 





























































































































Table: 4.19: Ni-40% Co annealed at 800°C for 15 minutes. 



Serial No. 

Component 

ID 

] 

Euler Angles 

fl[g) 

Volume % 


mUynHlH 


1 

{279}<7 11 
7> 

49.9 

37.4 

13.2 

4.61 

52.91 

2 

Cube 

0 

0.5 

0 

3.78 

7.69 

3 

Goss 

0 

45 

0 

2.5 

3.94 

4 

Retained S 

60.3 

35.3 

60 

2.1 

7.62 

5 

Bs Rex. 

70.3 

33.3 

50 

2.3 

16.65 


Table: 4.21: Ni-40% Co annealed at 800®C for 3 hr. 


Serial No. 

Component 

ID 

1 

Euler Angles 

m 

Volume % 

<i>i 


<1>2 

1 

Retained S 

55.5 

36.8 

61.1 

2.51 

24.3 

2 

Goss 

0 

45 

0 

2.05 

2.38 

3 

Cube 

0 

0.6 

0 

4.58 

15.84 

4 

{279}<711 

7> 

49.9 

39.2 

13.2 

4.25 

43.72 

5 

Bs Rex. 

69.4 

34.4 

50 

2.1 

12.31 




































































































Table: 4.22: Ni-40% Co annealed at 800°C for 10 hr. 


Serial No. 

Component 

ID 

] 

Euler Angles 

f(g) 

Volume % 


imilllRSHU 

<l>2 

1 

Goss 

4.8 

45 

0 

2.3 

5.9 

2 

Cube 

0 

0.6 

0 

3.3 

9.95 

3 

{279} <7 11 
7> 

51.7 

39.1 

12.4 

4.12 

51.07 

4 

Retained S 

64.2 

35.8 

60.2 

2.34 

16.88 

5 

Bs Rex. 

^ 75.5 

33.2 

45 

2.18 

14.49 


Table: 4.23: Ni-40% Co annealed at SOO^C for 20 hr. 


Serial No. 

Component 

ID 

1 

Euler Angles 

f(g) 

Volume % 

|||[||||||||^||^^ 



1 

Goss 

0 

45 

0 

2.37 

4.33 

2 

Cube 

0 

0.5 

0 

3.93 

7.02 

3 

{279}<7 11 
7> 

51.6 

39.6 

12.3 

4.4 

53.2 

4 

Retained S 

62.9 

35.3 

57.4 

2.71 

20.57 

5 

Bs Rex. 

69.9 

34.2 

• 50 

2.6 

13.16 


Table: 4.24: Ni-40% Co annealed at 800°C for 50 hr. 


Serial No. 

Component 

ID 

] 

Euler Angles 


Volume % 

||H|||||[|||^^ 

llllllllll^llllllllll^ 


1 

Cube 

0 

0.7 

0 

5.62 

11 

2 

{279}<7 11 
7> 

47.2" 

40.4 

13.5 

■ai 

48.48 

3 

Retained S 

60.9 

38.2 

59 

in 

6.46 

4 

Bs Rex. 

70.7 

34,3 

45 

2.34 

23.03 

5 

Goss 

0 

45 

0 

2.6 

9.77 
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Table: 4.25: Ni-60% Co annealed at 800°C for 15 minutes. 


Serial No. 

Component 

ID 

1 

Euler Angles 

f(g) 


<l>i 


HESmilil 


Goss 

0 

45 

0 

1.94 

2.15 

2 

{0 4 11}<8 

11 4> 

56.1 

21.1 

0 

1.95 

32.02 

3 

Cu 

90 

28.6 

45 

4.63 

22.76 



90 

35.8 

0 

4.9 

i 

|H 

5 

Near to Bs 

30 

45 

0 

1.7 

5.42 


Table: 4.26: Ni-60% Co aimealed at 800°C for 1 hr. 


Serial No. 

Component 

ID 

1 

Euler Angles 

m 

Volume % 

njRiHm 

______ 


1 

Near to 
rotated Goss 

19.2 

45 

0.2 

7.41 

HI 

2 

Rotated 

Goss 

90 

39.9 

0 

13.51 

47.07 

3 


54.1 

38.8 

10 

2.2 

4.55 

4 

Bs Rex. T1 

33.6 

41.6 

75 

4.2 

2.07 

5 

Bs Rex. 

81.2 

33.7 

30 

4.7 

17.9 


Table: 4.27: Ni-60% Co annealed at 800°C for 3 hr. 


Serial No. 

Component 

ID 

'] 

Euler Angles 

m 

Volume % 


llllllllllll^^ 

h. 

1 

Nearly 
rotated Goss 

15.5 

45 

0.1 

7.39 

37.4 

2 

Rotated 

Goss 

90 

39.7 

0 

11.4 

43.07 

3 

■ftlllETiVRa 

49.4 


10 ^ 

2 

3.59 

4 

Bs Rex. T1 



75 ^ 

3.9 

2.69 

5 

Bs Rex. 

wmnsm 

34 


3.3 

13.2 
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Table: 4.28: Ni-60% Co annealed at 800°C for 10 hr. 


Serial No. 

Component 

ID 

1 

Euler Angles 

%) 

Volume % 



4>2 

1 

Nearly 
rotated Goss 

15.3 

45 

0 

10.02 

42.42 

2 

{011}<655> 

44.8 

35.6 

3.2 

2.28 

25.28 

3 

Near to Bs 
Rex. T1 

30 

41.6 

75 

■■ 

1.41 

4 

Rotated 

Goss 

90 

45 

0 

10.7 

30.93 


Table: 4.29: Ni-60% Co annealed at 800°C for 20 hr. 


Serial No. 

Component 

ID 

] 

Euler Angles 

m 

Volume % 




1 

Rotated 

Goss 

90 

40 

0 

9.25 

29.56 

2 

Nearly 
rotated Goss 

15.1 

45 

0 

9.63 

23.4 

3 

iMiiBdLjg 

49.7 

39.3 

8.7 

2.03 

7.79 

4 

Near to Bs 
Rex. T1 

34.2 

41.5 

75 

4.29 

24.99 

5 

Goss 

0 

40 

0 

5.5 

6.45 

6 

Bs Rex. 

85.3 

33.1 

30 

2.4 

7.93 


Table: 4.30: Ni-60% Co annealed at 800®C for 50 hr. 


Serial No. 

Component 

ID 

] 

Euler Angles 

f(g) 

Volume % 




1 

Nearly 
rotated Goss 

18.3 

45 

0 

11.97 

40.63 

2 


90 

45 

0 

13.15 

36.96 

3 

Goss 

0 

45 

0 

6.5 

6.92 


KiHtBtfg 

49 

40.2 

5 

2.2 

3.17 

5 

Near to Bs 
Rex. T1 

30.2 

40.7 

75 

5.2 

4.85 

6 

Bs Rex. 

84.6 

35.2 

25 

2.8 

7.53 


107 







































































































Annealing Time (in hrs.) 


Figure 4.48; Variation of intensity of cube component as a function of annealing time 

for different samples. 
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Figure 4.49: Variation of volume fraction of cube component as a function of 
annealing time for different samples. 
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uiaug. u. luui cuoys. oul aner lu hours of annealing time the volume fraction 
of cube component for all the samples becomes nearly constant. 

The variation of intensity and volume fraction of cube component as a function of Co 
content for different annealing times are shown in Figures 4.50 and 4.51 
respectively. It can be observed from the figure that the intensity and volume fraction 
of the cube component decreases as the Co content is increased beyond 20 wt%. 

4.3.2: VARIATION OF CUBE COMPONENT AND ITS TWINS WITH 
ANNEALING TIME IN Ni-10% Co TO Ni-30% Co ALLOY 
The volume fractions of cube component and its twins, as a function of annealing 
time for different alloys, are shown in Figure 4.52. It can be observed from the figure 
that cube and its twin components constitute major volume fraction of samples. The 
volume fraction of cube and its twin component in Ni-10% Co possess almost 80% of 
the total volume from initial stages of annealing. Although there is no general trend 
for different samples but, if observed closely, it can be revealed that at the initial 
stages of annealing there is higher variation in volume fractions. With increasing 
annealing time the volume fraction become steady. 

The volume fraction of cube component and its twins, as a function of Co content for 
different annealing times are shown in Figure 4.53. Generally the volume fraction of 
cube and its twin component decreases if Co content increases beyond 20 wt%. 

4.3.3: VARIATION OF BsT COMPONENT WITH ANNEALING TIME IN Ni- 
10% Co TO Ni-30% Co ALLOY 

The variation of intensity and volume fraction of the BsT component {437}<4 1 1 7> 
as a function of annealing time for different alloys is shown in Figures 4.54 and 4.55 
respectively. The intensity of this component is nearly constant with annealing time 
and is very low. But there is a wide variation of volume fraction of this component 
upto 20 hours of annealing time after which it attains a constant value for different 
samples. In case of Ni-30% Co this component constitute about more than 40% of 
volume after 20 hours of annealing time. But the intensity of this conqwnent was 
relatively small which indicates that this component is quite spread out. 

4.3.4: VARIATION OF MAJOR TEXTURE COMPONENTS WITH 
ANNEALING TIME IN Ni-40% Co ALLOY 

The major texture components of Ni-40% Co alloy are cube conqjonent, BsR 

conqwnent {236}<385>, retained S component {123}<634>, and the {279}<7 11 7> 

component Among these «e have discussed about the cube component earlier. The 
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%1/M 



Wt% Co 


Figure 4.51: Variation of volume fraction of cube component as a fimction of Co 
content and subjected to different annealing times. 
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Annealing Time (in hrs.) 


Figure 4.52: Variation of volume fraction of cube and its twin con^nent as a 
function of annealing time for different samples. 
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Figure 4.54; Variation of intensity of BsT component {437}<4 1 1 7> as a function of 

annealing time for different samples. 
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Figure 4.55: Variation of volume fiaction of BsT component {437}<4 11 7> as a 
function of annealing time for different samples. 




intensity and volume fraction of rest of the components as a function of annealing 
time is plotted in Figures 4.56 and 4.57 respectively. The intensities are constant and 
nearly same for all the components. In case of volume fraction the {279}<7 11 7> 
component is nearly about 50% after 10 hours of annealing time. The plot of rest of 
the two components is mirror image of each other. It suggests that they grow at the 
expense of each other. Also considering both the plots it can be said that all three 
components are highly spread out. 

4.3.5: VARIATION OF MAJOR TEXTURE COMPONENTS WITH 
ANNEALING TIME IN Ni-60% Co ALLOY 

There are two major texture components of Ni-60% Co alloy. They are rotated Goss 
component {011}<110> and a component near to it {011}<114>. The intensity and 
volume fraction of this two components as a function of annealing time is plotted in 
Figures 4.58 and 4.59 respectively. The intensities of these two components are very 
nearly the same for all annealing times. The variations of the volume fractions are 
also similar for the two components and together they constitute almost 80% of the 
total volume. Again the low intensities of these components indicate that they are 
highly spread out. 
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Figure 4.56: Variation of intensities of different texture components in Ni-40% Co 

alloy as a function of annealing time. 
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Figure 4.57; Variation of volume fractions of different texture components in Ni-40% 

Co alloy as a function of annealing time. 





Figure 4.58; Variation of intensities of different texture components in Ni-60% Co 

alloy as a function of annealing time. 





Figure 4.59: Variation of volume fractions of different texture components in Ni-60% 

Co alloy as a function of annealing time. 




CHAPTER 5 
DISCUSSION 


Pure Ni has moderately high stacking fault energy. With upto 68% of Co, the Ni-Co 
alloys form continuous series of solid solution. This is evident from our study that 
throughout the annealing process the microstructures of all the samples consist of a 
single phase only. All the samples showed recrystallisation even after the shortest 
period of annealing time. Immediately after recrystallisation, grain growth started on 
further annealing. 

The frequency distributions of the relative energy ratios reveal that the 71/73 ratios of 
higher frequencies lie within values of 0.6 to 0 . 8 . In the class of 0.9 to 1 , 71/73 ratio 
drastically decreases, becoming zero at 71/73 “ !• The ratio of 71/73 can only be 1 when 
the dihedral angles become equal. At equilibrium the dihedral angles of 
polycrystalline materials come closer to 120° and the corresponding grain boundary 
energies also approach values close to one another. There is a wide range of 
frequencies of 7I/73 for different annealing time. This indicates that grain growth after 
recrystallisation is clearly a dynamic process with some grain boundaries becoming 
more active than others at a particular time and then some other boundaries may take 
over at a later stage. The profiles of the 71/73 plots, after the longest annealing time, 
namely 50 hours, do not show maxima at 71/73 = 1.0 for all the five Ni-Co alloys. This 
indicates that equilibrium is not yet achieved, otherwise this ratio should have gone to 
1 . 0 . 

It is observed from Figures 4.29 and 4.33 that the frequencies of the 72/73 ratios vary 
within a rather narrow range, eventually assuming maxima at 72/73 “ in no 

case the frequency of 72/73 is greater than 0.6 when 72/73 ” 1 * This result, coupled with 
the fact that 71/73 value never reached 1.0 for any of the alloys, it can safely be 
inferred that we are still far from reaching the equilibrium when at triple points the 
grain boundaries will be 120° ^art fixim one another. 

The texture results have been given in a tabular form in Tables 4.1 to 4 . 30 . It is clear 
from these tables that there are a number of components present in the textures of all 
the annealed samples of the five alloys. Ibe cube component and its twin orientations 
occupy the predominant volume fractions in Ni- 10 % Co to Ni- 30 % Co alloy. In Ni- 
40 % Co alloy the important components are the cube component, the BsR component 
{ 236 }< 385 >, the { 279)<7 11 7 > component, and the retained S component 



[ 123 }< 634 >. In Ni-60% Co alloy the impoitant comfionents are the rotated Ooss 

component {01 1 )<l 10> and the component {0U)<1 14> which is near to the rotated 
Goss component. 

After considering all the relative grain boundary energy parameters and texture 
parameters it was found that the best correlation can be found out if we compare the 
variation of the frequencies of yl/yS ratio as well as the variation of volume fraction 
of texture components with annealing time for different alloys. 

In case of Ni-I0% Co alloy it is seen that upto an annealing time of 20 hours there is 
no regular pattern in the variation of frequencies of yl/y3 ratio Avith annsaling time. 
This indicates a dynamic situation where some specific grain boundaries become 
active at one instant, which could change after further annealing, leading to some 
other specific boundaries becoming more active than the previous one. This kind of 
behavior has been obtained for all the Ni-Co alloys. Beyond 20 hours of annealing 
time the frequencies of different classes of yl/y3 ratio do not show any large variation. 
Also the same thing is true if we consider the variation of volume fraction of the cube 
component and the variation of volume fraction of cube and its twin component as a 
function of annealing time. 

For Ni-20% Co alloy even after 50 hours of annealing equilibrium is not reached both 
in terms of variation of frequencies of yl/y3 ratio with time and variation of volume 
fraction of cube and its twin component with time and variation of volume fraction of 
BsT component with time. 

In case of Ni-30% Co the frequency of the 0.6 to 0.7 class of yl/y3 ratio increases 
drastically after 20 hours of annealing time. This can be correlated with the increase 
of volume fraction of the BsT component with time after 20 hours of annealing. 

In case of Ni-40% Co alloy after 20 hours of annealing time the class of 0.7 to 0.8 of 
yl/y3 ratio decreases rapidly. From the texture plot of that alloy it can be seen that the 
volume fiiaction of the retained S component also decreases after 20 hours of 
annealing time. On the other hand there is an increase in the frequency of the classes 
of 0.5 to 0.6 and 0.6 to 0.7 of yl/y3 ratio after 20 hours of annealing. Similarly there is 
an increase in tbe volnme fiaction of (he BsR eomponeat wMi time after 20 hours of 
annealing. 

In case of Ni-60% Co alloy after 20 hours of annealing time the 0.7 to 0.8 and 0.8 to 
0.9 classes ofyl/yS ratio increase after 20 hours of annealing time. If we cmrsider the 



texture results the rotated Goss eomponeut and the {0U}<114> eomponent also 
increase in volume fraetion after 20 hours of annealing time. There is a rapid decrease 
of the class of 0,6 to 0.7 of r l/yS ratio after 20 honra of annealing. It can be eortelated 
with the decrease in the volume fraction of the eomponent {16 4 I3}<194> (neat to 
Bs rex. first generation twin) after 20 hours of annealing. 

The microstructural and textural study of Ni-Co alloys indicates that these materials 
are somewhat not stable, texture-wise and grain boundary energy wise even after 50 
hours of annealing at 800 C. It would be worthwhile to try to correlate texture and 
microstructure in these alloys by annealing them at 800“C for much longer periods of 
time or at temperatures higher than 800°C. 





1 . The recrystalHsalion textures of Ni-10% Co to Ni-30% Co alloys are very similar. 
These are comprised of mainly the cube, its twin components, and BsT 
component. 

2. In Ni-40% Co alloy the important recrystallisation texture components are the 
BsR component, the {279}<7 11 7> component, the retained S component, and 
the cube component. 

3. The recrystallisation texture of Ni-60% Co alloy consists mainly of rotated Goss 
component and the component {01 1 }<1 14>. 

4. There is a wide range of frequencies of y1/y 3 for different annealing times. This 
indicates that grain growth after recrystallisation is a dynamic process. 

5. The peaks of frequency distribution of yl/y^ ratio are spread out in the region of 
yl/y3=0.6 to yl/y3=0.8 for different annealing times. yl/Y3 ratio never reached 1.0 
for any of the alloys. 

6. In no case the frequency of y2/y3 is greater than 0.6 when y2/Y3=l . 

7. From 5 and 6 it can be concluded that equilibrium has not been reached within the 

annealing period used. 

8. Beyond 20 hours of annealing time the frequencies of different classes of yl/y3 
ratio do not show any large variation for Ni-10% Co alloy. This is not true for rest 

of the alloys. 

9. The frequency distribution of the grain boundary energy ratios correspond well 
with the volume fraction of the texture components in all the five alloys after 20 

hours of annealing treatment. 
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Appendix A 


Computer Programs 

Al. Program to read image data and convert into coordinate 
System 


/* mod exgb.c*/ 

#include <stdio.h> 

#include <conio.h> 

#define MAXPIXEL 15000 /* size for the max points analysed*/ 

#define MAXSEG 1000 /* maximum pixels in a gb segment */ 

#define GBPIXEL 'O' /* value of gb pixel */ 

int getap(int value[][2], int nc, int n, int pQ, int m); 
int test_tp(int value[][2], int pQ); 

main() 

{ 

intp[3],q[3],gb[MAXSEG],value[MAXPIXEL][2]; 
int i,j.k,l,nc,count; 
char c; 


for(i=0 j=0,nc=0;(((c=getc(stdm))!=EOF) && (nc<MAXPIXEL));) 
{ /* reading pixel coordinates */ 

if(c=^n') { /* end of line */ 

cprintf("(%i,%i)\n",i,j); 
i=i-H-; 
j=0; 


} 

ifi:c = GBPIXEL) { 

value[nc][0] = j; /* array stores x,y position */ 

value[nc][l] = i; of gb pixel */ 

nc-H-; 


} 

if((c = '0') ||(c=='l'))j-H-; 

(nc>=MAXPIXEL) cprintfCWaming: Could not read all of image data!\n"); 


/* set triple point counter */ 


count = 0; 
for (i=0;i<nc;i++) 

{ 

k = getap(value,nc,i,p,3); ^ 

if (jt == 3) { /* 1 p“ 3 : it is a triple point / 

if(test_tp(value,p)=l) continue; 




/* end of line marker */ 


for (j=0;j<3u-H-) 

{ > 

1 = 2 ; 

gb[0] = i; 

gb[l]=p[j]; 
do { 

if ((k = getap{value,nc,gb[l-l],q,3)) = 2) { 
if (gb[l-2] == q[0]) /* exclude the duplicate pixel */ 

gb[l] = q[l]; 

else 

gb[l] = q[0]; 

1-H-; 

} 

} while ((k == 2) && (1 < MAXSEG)); 

fbr(k=Od!;<l;k-H-) printl("%i %i ",value[gb[k]][0],value[gb[k]][l]); 

if (j < 2) printf("- l\n"); /* end of line marker */ 

} 

count-H-; 

} 

} 

printi(" %d\n", count); 

printf("-99\n"); /* end of file marker */ 

} 


int getap(int value[][2], int nc, int n, int pQ, int m) 

{ /* returns no. of adjacent pbcels touchii^ pixels n */ 

/* on return, p[] contains the first m adjacent points */ 

int i,j,k,x,y,dx,dy; 


X = value[n][0]; 
y = value[n][l]; 
i = n; 

j = -i; 

k = 0; 
do{ 
i = i+j; 

if (i >= nc) break; 
if(i<0){ 
i = n; 

j = i; 

} 

else { 

dy = value[i][l] - y; 
if((dy<-l)|l(dy>l)){ 

ifO = -l){ 

i = n; 

j = i; 

} 

else 


/* no more' gb pixels to check */ 


/* rMWMi tn fiirtha’ */ 


} 

else { 

dx = value[i][0] - x; 
if ((dx > -2) && (dx < 2)) { 

if (k < m) p[k] = i; /* found a pixel "i" touching pixel "n" */ 
k-H-; 


} 

} 

} 

} while(l); 
retum(k); 

} 


int test_tp(int value[][2], int pQ) 

{ 

int x[3],y[3]; 
int i,j; 

for(i=0;i<3;i-H-) 

{ 

x[i]=value[p[i]][0]; 

y[i]=value[p[i]][l]; 

} 

for(i=0;i<3;i-i-i-) 

{ 

j=i+i; 

ifG=3)j=0; 

if(x[i]=xD]) 

{ ifl:y[i]==(yD]+i)lly[il=(yD]-i)) 

retum(l); 

} 

if(y[i]=y[j]) 

{ if(x[i]=(xD]+l)||x[i]=(xQ>l)) 
retum(l); 

} 

} 

retum(0); 

} 



A.2 Program for Evaluaring Plane Dihedral Angles 

/* modtpaB.c*/ 

#include <ermo.h> 

# inc lude<std io . h> 

#include<alloc.h> 

#include<math.h> 

#include<conio.h> 

#include "modcfitS.c" /* this file is given in Appendix A.3 */ 

#define MAX_DATA 100 
Mefine MAX_TP 400 

#define DISTANCE(xl,yl,x2,y2) ((xl-x2)*(xl-x2)+(yl-y2)*(yl-y2)) 

typedef struct { 

int order,npts; 
float x,y,r2,*a; 

} polycoeff; 

typedef struct { 

float tx,ty,beta[3]; 
polycoeff g[3]; 

} triplepoint; 

int tpfit( triplepoint tQ, int *ntp); 

int gb_line( float x[], float y[], int *ndat); 

void pda(triplepoint tQ, int ntp); 

get_distance(double d[], double si, double s2, triplepoint t, int i); 
double slope( float aQ, int order, int x); 

mainQ 

{ 


int ntp; 

triplepoint t[MAX_TP]; 
cprintf("\nProgram Begin\n"); 
tpfit(t,&ntp); 
pda(t,ntp); 

cprintf("ntp = %i\n",ntp); 


void pda(triplepoint tQ.int ntp) 

{ 

int k,ij,p,q,r; 
double s[3],d[3]; 

FILE 



if((fp ~ fopen( coeff.log , w )) NULL) cprintf("unable to create coefF.logW); 
fprintf(fp,"#,tpx,tpy,order,aO,al ,a2,a3,a4,a5,a6\n"); 

printf("#,tpx,tpy,ordl,r2_l,ord2,r2_2,ord3,r2_3,betal,beta2,beta3,sumbeta\n"); 

for(k=0;k<ntp;k-H-) 

{ 

if(t[k].g[01.npts<=4|it(k].g[l].npts<=4|lt[k].g[2].npts<=4) 

{ 

/* printf("\n%f %f has too few points to continue\n",t[k].tx,t[k],ty);*/ 
continue; 

} 

if(t[k].tx=-l) 

continue; 

printf("%i,%f,%f,",k+l,t[k].tx,t[k].ty); 


for(i=0;i<3;i-H-) 

{ 

s[i]=slope(t[k].g[i].a,t[k].g[i].order,t[k].tx); 
/* printf(" slope %d is %lf ",i,s[i]);*/ 

} 

I* printf("\n");*/ 


for(i=0;i<3;i++) 

{ 

j=i+l; 

ifa=3)j=0; 

get_distance(d,s[i],s[j],t[kl,0; 

t[k].beta[i]=acos((d[0]+d[l]-d[2])/(2*sqrt(d[0])*sqrt(d[l])))*(l 80/3. 1412), 

} 

check = ' '; 

if((t[k].beta[ 0 ]+t[k].beta[l ]+t[k].beta[2])<360) { 
check = 

for(p=0;p<3;p++){ 


q=p+l; 
if(q=3) q=0; 
r=q+l; 

if(((t[k].beta[p]+t[k].beta[q3)>=KtEk].bete 

&&((t[k].beta[p]+t[k].beta[q])<=(t[k].beta[r]+0.001))) 


t[k].beta[p]= 180 -t[kl.beta[p]; 
t[k].beta[q]=l 80-t[k].beta[q]; 


} 


} 


} 



for(p=0;p<3;p++) prmttr’%i.%f,",t[k].g[p].order,t[k].g[p].r2); 
for(p=0;p<3 ;p-H-) printfi:"%f,",t[k].beta[p]); 
printf("%r,t[k].beta[0]+t[k].beta[l]+t[k].beta[2]); 
if (check = ' ') 
printf("\n"); 

else 

printfl;",%c\n", check); 

for(p=0;p<3;pHH-) 

{ 

fprintf(fp,"%i,%f,%f,%i",k+ 1 ,t[k].tx,t[k].ty,t[k].g[p].order); 
for(i=0;i<=t(k].g[p].order;i-H-) fprintf(fp,",%f',t[k].g[p].a[i]); 
fprintf(fp,"\n"); 

} 

} 

fclose(fp); 

} 

get_distance( double d[], double si, double s2, triplepoint t, int i) 

{ 

double c 1 ,c2,y 1 ,y2,xl ,x2,ml ,m2,kl ,k2; 
intj; 

j=i+i; 

if(j=3)j=0; 

cl==t.ty-sl*t.tx; 

c2=t.ty-s2*t.tx; 

xl=t.g[i].x; 

x2=t.g[j].x; 

yl=t.g[i].y; 

y2=t.gD].y; 

if (si = 0) { 
xl = t.g[i].x; 
yl = t.ty; 

} - 
else { 

ml=-(l/sl); 

kl=(yl-(ml*xl)); 

xl=((kl-cl)/(sl-ml)); 

yl=(((ml*cl)-(sl*kl))/(ml-sl)); 

} 

if(s2 = 0) { 
x2 = t,g0].x; 
y2 = t.ty; 

} 

else { 

m2=-(l/s2); 

k2=(y2-(in2*x2)); 



y2=(((ni2*c2)-(s2*k2))/(ni2-s2)); 

} 

d[0]=DISTANCE(t.tx,t.ty,x 1 ,y 1 ); 
d[ 1 ]=DISTANCE(t.tx,t.ty,x2,y2); 
d[2]=DISTANCE(xl,yl,x2,y2); 

} 


double slope( float a[], int order, int x) 

{ 

int i; 

double sum = 0; 

for(i= 1 ;i<=order;i+4-) 
sum=sum+i*((double) a[i])*pow(x,i-l); 

retum(sum); 


} 

int tpfit(triplepoint t[], int *ntp) 

{ 

float x[MAX_DATA],y[MAX_DATA],sig[MAX_DATA],a[MAX_SIZE]; 
int ndat,i,ordj,l,err,fe,upperord; 
float r2; 

/* FILE ^fp;*/ 

/* DATA INPUT & CALCULATIONS */ 

for (i=0;i<MAX_DATA;i-H-) sig[i] = 1; 

/* fp = (FILE *) fopen(feame,"r"); 
if(fp = NULL){ 

printfi;''file error %iAn'',ermo); 
exit(O); 

} 

*1 

(*ntp) = 0; 
do 
{ 

for(i=0;i<3;i++) 

^ fe = gb Iine(x,y,&ndat); 

err = polyfit(x,y,sig,ndat,l,6,&ord,a,&r2,l,x[0],y[0]); 

if(!err) 

t[*ntp].tx = x[0J; 
t[*ntp].ty = yE01; 


t[*ntp].g[i].y=y[ndat-l]; 
t[*ntp].g[i].r2 = r2; 
t[*ntp].g[i].npts=ndat; 
if ((t[*ntp].g[i].a = (float*) 

^ malloc((ord+ 1 )*si 2 eof(float))) = NULL) 

cprintf(''waming: out of memory in tpfit\n"): 
getchO; 

^ retum(-l); /* out of memory */ 

^ for(j=0d<=ord;j++) t[*ntp].g[i].a[j] = aQ]; 
else { 

t[*ntp].tx = -l; 
for(l=i+l;l<3;I++) 

fe = gb_line(x,y,&ndat); 

break; 

} 

} 

(*ntp)-H-; 

} while(fe != -99); 
retum(0); 

} 

int gb_line( float x[], float y[],int *ndat) 

{ 

int num; 


(*ndat) = 0; 
do 
{ 

scanf("%i",&num); 
if(num >= 0) { 

if(*ndat<MAX_DATA) 

x[*ndat] = (float) num; 
scanfl["%i",&num); 
if (♦ndat < MAX_DATA) { 

y[*ndat] = (float) num; 
(*ndat)++; 

} 

} 

} while(num >= 0); 
retum(num); 


} 



A.3 Program for Curve Fitting 

/* modcfit.c*/ 

#include <stdio.h> 

^include <stdlib.h> 

^include <math.h> 

^include " cfit.h" /* This header file is given in Appendix A.4*/ 
/* 

main() 

( 


float x[100], y[100], sig[100], a[MAX_SIZE], xi; 
int i, ord; 


for (i=0;i<50;i-H-) 

{ 

xi = ((float) i)/50.0; 
x[i] = xi; 

y[i] = 0.01+0.02*xi+0.04*xi*xi + ((float) random(100))/100; 
sig[i] = 1; 

} 

polyfit(x,y,sig,50, 1 ,7,&ord,a,0,0,0); 

} 

*/ 

int polyfit( float x[],float y[], float sig[],int ndat,int ordl,int ord2,mt *ord,float a[],float 
*r2,int ianch,float anch_x,float anch_y) 

{ 

/* Polynomial Least Square Fit: 

X, y: "ndat" data points; sig: "ndat” standard deviations; 
ordl,ord2: search poly order ordl to ord2 and return best value in *ord (to 
implement) 

a: on return contains "*ord+l" polynomial coefficients 
ianch: non-zero value => set anchor at "pch_x,anch_y" 
return value: non-zero => no solution possible */ 

int temp,ij,degl,deg2,err, ia[MAX_SIZE]; 
float chisq.chisq 1 ,al [MAX_SIZE],ftest; 
double yc,ym; 

long double covar[MAX_SIZE][MAX_SIZE], a_x, st, sr; 
void (♦fh)(float, long double Q, hit); /* function pointer •/ 
if (ord2<ordl) SWAP(ordl,ord2) 

if ((ordl<0)l|(otd2>MAX_ORDER)) retutn(-l); /♦ invalid polynomial order V 
if (ndat<==ordl) retum(-l); /* invalid no. of data points */ 
if (ndat<=ord2) ord2=ndat-l; /* reduce the upper polynomial order / 



fh = polyfuncl; 

if (ianch) { /* set anchor point *! 

anchor_x = anch_x; 
anchor_y = anch_j; 
fh = polyfunc2; 

} 

for( *ord=ord 1 ;*ord<=K)rd2;(*ord)-H-) 

{ 

for(i=0;i<=*ord;i-H-) ia[i] = 1; 
if (ianch) ( /* set anchor for each order value */ 

ia[0] = 0; 
a[0] = anchor_y; 

} 

err = gllsq(x,y,sig,ndat,a,ia,(*ord)+l,covar,&chisq,fh); 
if (err != 0) retum(err); /* Error: solution not possible */ 
if (ianch) ( 
a_x = 1 ; 

for (i=l ;i<=*ord;i-H-) 

{ 

a_x = a_x*((long double) anchor_x); 
a[0] -= a[i]*a_x; 

} 

} 

if((*ord)>ordl) 

{ 

degl=ndat-(*ord); 
deg2=ndat-((*ord)+ 1 ); 
if (chisql < chisq) break; 
if((deg2!=0)&&(chisql !=chisq)) 

flest=(((deg 1 -deg2)*chisq)/((chisq 1 -chisq)*deg2)); 
else 
break; 
if(ftest >1) 
break; 

} 

for(i=0;i<=(*ord);i-H-) 
al[i]=a[i]; 
chisq l=chisq; 

} 

(*ord)=(*ord)-l; 

for(i=0;i<=(*ord);i++) a[i] = al [i]; 

ym = 0; /* calculation of correlation coefficient */ 

for(i=0;i<ndat;i-H-) ym += y[i]; 

ym = ym/ndat; 

sr = St = 0; 

for(i=0;i<ndat;i++) 

^st += (long double) ((y[i] - ym)*(y[i] - ym)); 


> 


for(i-Oa<=*ordJ-H.) yc += aUl*pow((double) >c[l), (double) i). 
sr += (long double) ((y[i] - yc)*(y[i] - yc)); ^ ’ 


if((st = 0.0) II (sr = 0)) 
(*r2)=l; 


else 

(*r2) = (float) ((st - sr)/st); 
retum(O); 


} 


void polyfuncl (float x, long double afuncQ, int ma) 

/* returns basis functions (at x) in "aflinc" in a polynomial of order "ma-l" */ 
int i; 

alunc[0] = 1; 
for (i=l;i<ma;i-H-) 

{ 

afuncfi] = ((long double) x) * afunc[i-l]; 

} 

} 


void poIyfiinc2(float x, long double afiincQ, int ma) 

{ 

/* returns basis functions (at x) in "afiinc" in a polynomial of order "ma-l" 
under constraint that curve is anchored at: (anchor_x, anchor_y) */ 

int i; 

long double xl = 1, x2 = 1; 


afunc[0] = 1; 
for (i=l ;i<ma;i-H-) 

{ 

xl =xl*x; 

x2 = x2*((long double) anchor_x); 
afunc[J] = xl - x2; 

} 

} 

int gllsq(float x[], float y[], flpat sig[], int ndat, float a[], int iaO, int ma, long double 
covar[][MAX_SIZE], float *chisq, void (♦gfunc)(float, long double Q. int)) 

{ 

/* General Linear Least Square Routine: 

X, y: "ndat" data points; sig: "ndat" standard deviations; ma: # of parameters; 
ia: "ma" size with only zero or non-zero values; zero talue implies, corre^nding 
parameter held fixed to value stored in "a"; gfunc: user routine to su|^ly "ma" 
basis functions evalated in array afiinc (size: "ma"); 
returns ChiSquare in "chisq" and covariance matrix in covar; 
return value: non-zero — > no solution possible */ 



int err,i,j,k,l,m,mfit=0; 
float ym,wt,sum,sig2i; 

long double beta[MAX_SIZE],afunc[MAX_SIZE]; 

for (j=0:j<may-M-) 
if (iaUl) mfit-H-; 

if (mfit == 0) retum(- 1 ); /* no parameters to be fitted */ 


for (j=0;j<mfity-H-) 

{ /* Initialize the (symmetric) matrix */ 

for (k=0;k<tnfitdc-H-) covar[j][k]=0.0; 
beta[j]=0.0; 

} 

for (i=0;i<ndat;i-i-+) 

{ /* Loop over data to accumulate coeffcients of the normal equations */ 

(*gfimc)(x[i],afunc,ma); 
ym=y[i]; 
if (mfit < ma) 

{ /* Subtract off dependences on known pieces of the fitting function */ 

for (j^0y<may++) 
if (liafj]) ym -= a[j]*afunc[j]; 

} 

sig2i=l.0/SQR(sig[i]); 
for (j=-l,l=0;l<ma;l-H-) 

{ 

if (ia[l]) { 
wt=afunc[l]*sig2i; 

for (j-H-,k=-l,m=0;m<=l;m-H-) 
if (ia[m]) covar[j][-H-k] += wt*afunc[m]; 
beta[j] += ym*wt; 

} 

} 

} 

for 0=1 j<mfitu++) /* Fill in above the diagonal from symmetry */ 
for (k=0;k<j^c-H-) 
covar[k]0]=covar[j][k]; 

if ((err=gaussj(covar, mfit, beta)) != 0) /* Matrix solution */ 

retum(err); /* no solution possible */ 
for 0=-l ,l=0;l<ma;l++) 

if (ia[l]) a[l]=beta[-H-j]; /* Partition solution to appropriate coefifecients a */ 
*chisq=0.0; /* and Evaluate ChiSquare of the fit */ 

for (i=0;i<ndat;i-H') 

{ 

(*gfunc)(x[i],afimc,ma); 
for (sun^.0j=0y<may-H-) 
sum += a[j]*afunc[n5 
♦chisq += SQR<(y[i]-sum)/sig[i]); 

Lvsrt(covarwia,iifit); /• Sort covariaiKK imtrix to true order of King cowfeiits 



retiim(O); 

} 

void covsrtdong double covar[][MAX_SIZE], int ma, int ia(], int mfit) 

/* Expand in storage the covariance matrix covar, so as to take into account 
parameters that are being held fixed (For latter, return zero covariances) */ 

int i,j,k; 

long double temp; 
for (i=mfit;i<ma;i++) 

for (j=Oo<=iy++) covar[i][j]=covarG][i]=0.0; 
k=mfit; 

for (j=ma-ly>=0u-) { 

if (iaUl) { 

for (i=0;i<ma;i++) SWAP(covar[i][k],covar[i][j]) 
for (i=0;i<ma;i-H-) SWAP(covar[k][i],covar[j][i]) 
k-; 

} 

} 

} 

int gaussj(long double a[][MAX_SIZE], int n, long double b[]) 

{ 


/* Gauss-Jordon Method; uses Full Pivoting 
a: square matrix of size "n x n”; b: vector of size n 
on return: "a" contans inverse of "a" and the solution in "b" 
return value: == 0 => no error; !== 0 => no solution possible */ 

intindxc[MAX_SIZE],indxr[MAX_SIZE],ipiv[MAX_SIZE]; 

int i,icol,irowj,k,l,U; 

long double big,dum,pivinv,temp,templ; 

for (j=0y<ny++) ipiv(J]=0; /* set pivot tracking elements to zero */ 

for (i=0;i<n;i++) /* main loop for pivot search */ 

{ 

big=0.0; 

for (j=0y<nu++) 

^ if (ipivjj] != 1 ) /* skip, if di^onal element in jth row is a pivot */ 

{ 

for (k= 0 dc<nd£-H-) 

{ 

if (ipiv[k] == 0) 

{ 

if(an][kl<0) 



else 


tempi = a|j][k]; 
if (tempi >=big) 

{ 

big = tempi; 
irow=j; 
icol=k; 


else if (ipiv[k] > 1) retum(l); /* Error: Singular Matrix-1 */ 

} 

} 

} . 

■H-(ipiv[icol]); /* irow, iool should contain the position of pivot */ 

if (irow != icol) /* if irow = icol, then pivot already in place */ 

{ /* else position pivot by swapping rows */ 

for (l=0;l<n;l-H-) SWAP(a[irow][l],a[icol][l]) 

SWAP(b[irow],b[icol]) 

} 

indxr[i]=irow; /* keep track of the original pivot position */ 

indxc[i]=icol; 

if (a[ icol] [icol] — 0.0) retum(2); /* Error: Singular Matrix-2 */ 

pivinv=l .0/a[icol][icol]; /* Now divide the pivot row by the pivot element */ 
a[icol][icol]=L0; 

for (l=0;l<n;l-H-) a[icol][l] *= pivinv; 
b[icol] *= pivinv; 


for (II=0;Il<n;ll++) /* Reduce the rows except the pivot one *1 

{ 

if (11 != icol) 

{ 

dum=a[ll][icol]; 

a[ll][icol]=0.0; 

for (l=0;l<n;l-H-) a[Il][l] a[icol][l]*dum; 
b[U]-=b[icol]*dum; 

} 

} /* End of Main Loop *1 

for (l=n-l ;1>=0;1“) /* Unscramble the solution by interchanging columns ♦/ 

{ 

if(mdxr[l] !=indxc[l]) 
for (k=0;k<n;k-H-) 

SWAP(a[k][indxr[l]],a[k][indxc[l]]); 

} 

retum(0); 

} 



A.4 Header File Used by Curve Fitting Program 

/* cfit.h*/ 

#define SWAP(a,b) {temp=(a);(a)=(b);(b)=temp;} 
static float sqrarg; 

#define SQR(a) ((sqrarg=(a)) = 0.0 ? 0.0 : sqrarg*sqrarg) 

#define MAX_ORDER 9 

#define MAX_SIZE MAX_ORDER+l 

static float anchor_x, anchor_y; 

int polyfit(float xfj.float y[],float sig[],int ndat,int ordl,int ord2,int *ord,float a[], float 

*r2,int ianch,float anch_x,float anch_y); 

void polyfimcl (float x, long double afuncj], int ma); 

void polyfunc2(float x, long double afunc[], int ma); 

void covsrt(long double covar[][MAX_SIZE], int ma, int ia[], int mfit); 

int gllsq( float x[], float y[], float sig[], int ndat, float a[], int ia[], int ma, long double 

covar[][MAX_SIZE], float chisq[], void (*gfunc)(float, long double [], int)); 

int gaussj(long double a[][MAX_SIZE],int n, long double b[]); 



A.5 Program to Convert PDA to TDA 

/*pda2lda.c*/ 

#include<stdio.h> 

#include<stdlib.h> 

#include " gaussj.c" /* This program is give in appendix A.6*/ 

main() 

{ 


FILE *fp,*fa,*fb; 

long dpuble a[MAX_SIZE][MAX_SIZE]; 
long double b[MAX_SIZE]; 
char fhamel[50],fiiame2[50]; 
int i=0,j=0,nc; 
float x; 

fp=fopen("table.txt","r’'); /* Table 2.3 is used here*/ 

if(fp=NULL) 

printf("errorl"); 

do { 

fscanf( fp,"%f ',&x); 
iffx < 0) 

( 

i=i+l; 

j=0; 

} 

else { 

a[i][j] = (long double) x; 

j=j+i; 

} 

} while(x > -99); 

nc = i; 
fclose(fp); 

printf("Input File Name: "); 
scanf("%s",fiiame 1 ); 
print:^"Output File Name: "); 
scanl^"%s'',fiiame2); 
fe=fopen(fiiame 1 ,"r"); 
fb=fopen(fiiame2,"w"); 
ij9:fa=NULL) 

{ 

printfl^"error2\n''); 

exit(O); 



for (i=0;i<nc;i-H-) 

{ 

fscanti:fa,"%Lr,&bti]); 

} 

fclose(fp); 

fclose(fa); 


gaussj(a, 1 7,b); 

for(i=0;i<17;i-H-) 

{ 

printf("%LAn",b[i]); 

fprintf(fb,"%Lf\n",b[i]); 

} 

getchO; 

} 



A.6 Program Used by PDA to TDA Conversion Program 

/*gaussj.c*/ 

#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include cfit.h /* This header file is given in appendix A.4*/ 

/*main() 

{ 


float x[100], y[100], sig[100], a[MAX_SIZE], xi; 
int i, ord; 


for (i=0;i<50;i-H-) 

{ 

xi = ((float) 0/50.0; 
x[i] = xi; 

yfi] = 0.01+0.02*xi+0.04*xi*xi + ((float) random(100))/100; 
sig[i] = l; 

}*/ 

polyfit(x,y,sig,50, 1 ,7,&ord,a,0,0,0); 


} 

int gaussj(long double a[][MAX_SIZE], int n, long double b[]) 

{ 


/* Gauss- Jordon Method; uses Full Pivoting 
a: square matrix of size "n x n”; b: vector of size n 
on return: "a" contans inverse of "a" and the solution in "b” 
return value: = 0 => no error ; !— 0 => no solution possible */ 
int indxc[MAX_SIZE],indxr[MAX_SIZE],ipiv[MAX_SIZE]; 

int i,icol,irowj4c,I,U; 

long double big,duni,pivmv,temp,templ; 

for (j=0;j<n;j-H-) ipiv0]=0; /* set pivot tracking elements to zero */ 


for (i=0;i<n;i-H-) 

{ 

big==0.0; 

for 0=0y<ny++) 

{ 

if(ipiv0] != 1) 

{ 

for (k=0;k<n;k-H-) 

{ 

if(ipiv[k) = 0) 


/* main loop for pivot search */ 


/* skip, if diagonal element in jth row is a pivot */ 



if (aO][k] < 0) 

tempi =-a[j][k]; 

else 

tempi = a[j][k]; 
if (tempi >= big) 

{ 

big = temp 1 ; 
irow=j; 
icol=k; 

} 

else if (ipiv[k] > 1) retum(l); /* Error: Singular Matrix-1 */ 

} 

} 

-H-(ipiv[icol]); /* irow, icol should contain the position of pivot */ 

if (irow != icol) /* if irow = icol, then pivot already in place */ 

{ /* else position pivot by swapping rows */ 

for (l=0;l<n;l++) SWAP(a[irow][l],a[icol][l]) 

SWAP(b[irow],b[icol]) 

} 

indxr[i]=irow; /* keep track of the original pivot position */ 

indxc[i]=icol; 

if (a[icol][icol] = 0.0) retum(2); /* Error: Singular Matrix-2 */ 

pivinv= 1 .0/a[icol] [icol] ; /* Now divide the pivot row by the pivot element */ 

a[icol][icol]=1.0; 

for (l=0;l<n;l-H-) a[icol][l] *= pivinv; 
b[icol] *= pivinv; 


for (ll=0;ll<n;lI-H-) /* Reduce the rows except the pivot one */ 

{ 

if (11 != icol) 

{ 

dum=a[ll][icol]; 

a[ll][icol]=0.0; 

for (l=0;l<n;l-H-) a[ll](l] -= a[icol][l]*dum; 
b[ll] -= b[icol]*dum; 

} 

I ^ /* End of Main Loop */ 

for (l=n- 1 ;1>=0;1--) /* Unscramble the solution by interchanging columns ♦/ 

{ 

if (indxr[l] != indxc[l]) 
for (k=0dt<n;k++) 

SWAP(a[k][indxr[l]],a[k][indxc[I]]); 

} 

retum(O); 



A.7 Program to Calculate Grain Boundary Energy Ratios 

/♦gdist.c*/ 

#include<stdio.h> 

#include<math.h> 

#include<stdlib.h> 

#define PI 3.1415927 

float gdist( float class[], float freq[], float cum[],int n); 

main() 

{ 

FILE *fe,*fb,*fc,*fd,*fe; 

int ij,m; 

float 

class[50],freq_max[50],freq_min[50],cuml [50],cum2[50],al ,a2.a3,bl ,b2,b3,g_max.g 

_min,g_middle,g_min_max,g_tniddle_max; 

char max_tda[50],min_tda[50],outputl [50],output2[50]; 

fa=fopen("\\tc\\source\\arijit\\c lass.txt", "r"); 

for(i70;i<17;i-H-) 

{ 

fscanf(fa,"%f’,&class[i]); 

} 

printf("\nlnput MAX_TDA File : "); 
scanf('’%s",max_tda); 
printfCVnlnput MIN_TDA FUe : "); 
scanf(''%s",imn_tda); 
print^"\nOutput File (gl/g3):"); 
scanfl[”%s", output 1 ); 
print^"\nOutput File (g2/g3):"); 
scanfl^"%s",output2); 

fb=fopen(max_tda,"r"); 
fc=fopen(min_tda,"r"); 
fd=fopen(output 1 ,"w"); 
fe=fopen(output2, " w”); 

for(i=0;i<l7;i-H-) 

fscaniR:fb,"%f’,&freq_max[i]); 

for(i=0;i<17;i'H-) 

fscan^fc,"%f',&freq_min[i]); 


printf("No of Iteration : "); 


forG==- Ou<m;j++) 

{_ 

al -gc3ist(class,freq_max,cuml , 1 7); 

a2=gc:i ist(class,freq_min,cuni2, 1 7); 

a3=3<^0-(al+a2); 

bl=si^(PI*al/180); 

b2=si^(PI*a2/180); 

b3=si^(PI*a3/180); 

iflbl=-b2) 

{ 

if(b2:==>b3) 

{g_»rnax=bl ;g_min=b3;g_middle=b2;} 
else 

if(b» l>b3) 

{g_max=bl ;g_min=b2;g_middle=b3;} 

els^s 

{g_max=b3;g_min=b2;g_middle=bl ;} 

} 

else 

{ 

if(b3:::>b2) 

{ g max=b3 ;g_min=b 1 ;g_niiddle=b2;} 

else 

if(t=>l<b3) 

{g_max=b2;g_min=bl ;g_middle=b3;} 

eis^ 

{g_raax=b2;g_inin=b3;g_middle=bl ;} 

} 

g_mii:m-_max=g_inin/g_max; 

g_inicl cile_max=g_middle/g_max; 

if ((g min_max < 0) || (g_imddle_max<0)) 

{ 

print^”%f, %£, %C %f\n",al,a2,g_nim_max,g_middle_max); 

} 

fprintfTfd,''%f\n",g_min_max); 

} 

fcIoseCfa)^ 

fcloseCfb); 

fcloseCfc); 

fcloseCf^)J 

getch(7>; 

} 

float g^ist(float classO, float freq[],float cumn,int n) 

{ 



/* calculation of cumulative frequency*/ 


ciim[0]=freq[0]; 
for(i=l ;i<n;i-H-) 

{ 

cum[i]=cum[i-l ]+freq[i]; 

} 

f=((float)rand())/RAND_MAX; /* random number generation*/ 
i=0; 

while(f > cum[i]) 
i++; 

upper=class[i]; 

uf=cum[i]; 

if(i>0) 

{ 

lower=class[i-l]; 

ljN:um[i-l]; 

} 

else 

{ 

lower=0; 

lf=0; 

} 

if(uf>lf) 

{ 

y=(upper* (f-If)+lo wer*(uf-f))/(uf-lf) ; 
retum(y); 

} 

else 

retum(-l); 

} 
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